Fast Hamiltonicity checking via bases of perfect matchings 
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Abstract 



For an even integer t > 2, the Matching Connectivity matrix H t is a matrix that has rows 
and columns both labeled by all perfect matchings of the complete graph K t on t vertices; 
an entry Ht\Mx,M^[ is 1 if Mi U Mi is a Hamiltonian cycle and otherwise. Motivated 
by the computational study of the Hamiltonicity problem, we present three results on the 
structure of T-L t : We first show that % t has rank at most 2*/ 2_1 over GF(2) via an appropriate 
factorization that explicitly provides families of matchings X f forming bases for %t . Second, 
we show how to quickly change representation between such bases. Third, we notice that 
the sets of matchings X t induce permutation matrices within H t . 

Subsequently, we use the factorization to obtain an 1.88S""nP^ time Monte Carlo al- 
gorithm that solves the Hamiltonicity problem in directed bipartite graphs. Our algorithm 
as well counts the number of Hamiltonian cycles modulo two in directed bipartite or undi- 
rected graphs in the same time bound. Moreover, we use the fast basis change algorithm 
from the second result to present a Monte Carlo algorithm that given an undirected graph 
on n vertices along with a path decomposition of width at most pw decides Hamiltonicity 
in (2 + v / 2) pw « c ' < ' 1 ' time. Finally, we use the third result to show that for every e > this 
cannot be improved to (2 + y/2 — e)^n ^ time unless the Strong Exponential Time Hy- 
pothesis fails, i.e., a faster algorithm for this problem would imply the breakthrough result 
of a (2 - e) n time algorithm for CNF-Sat. 

1 Introduction 

The Hamiltonicity problem and its generalization to the traveling salesman problem are widely 
acknowledged to be two of the most famous NP-complete problems. Many classical algorithms 
were invented to tackle these problems (and variants thereof), among them Cristofides' approx- 
imation algorithm [12], the Lin-Kernighan heuristic [31] . and a polynomial time approximation 
scheme for Euclidean TSP [T]. 

A very early and classical result belonging to this list is due to Bellman [U [5] and, indepen- 
dently, Held and Karp [23]; it demonstrates that the traveling salesman problem can be solved 
in 0(n 2 2 n ) time, where n denotes the number of vertices of the input graph. In order to get this 
result they introduced dynamic programming over subsets, which became a fundamental design 
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paradigm for obtaining exact exponential time algorithms (see for example |17j and \20\ Chap- 
ter 3]). Considerable effort has been taken to improve these algorithms: More space efficient 
algorithms were already given independently by several authors [31 [271 130] . 

In his influential survey, Woeginger [38] brought renewed attention to a question suggesting 
itself already from the 60's: Can either the traveling salesman problem or the Hamiltonicity 
problem be solved in (2 — e) n n°^ for some e > 0? An affirmative answer was given for the 
special cases of bounded degree [181 UJ [261 [21] an d claw-free graphs [11]. In a breakthrough result, 
Bjorklund [6] partially resolved the open question by giving an l.QQ n rP^ time Monte Carlo 
algorithm deciding undirected Hamiltonicity, using as a first step time algorithm 

for bipartite graphs. Unfortunately, it seems hard to derandomize this result or to extend it 
to the traveling salesman problem without incurring a pseudo-polynomial dependence on the 
input weights. The algorithm of [6J instead of keeping track of all the vertices along the so 
far constructed path (as Bellman, Held and Karp), labels only some vertices and some edges 
used on the path. Bjorklund's algorithm uses counting modulo two as a tool for canceling 
unwanted self-crossing walks, however it relies on the assumption that the Hamiltonian cycles 
under consideration satisfy particular counting invariants concerning the number of vertices and 
edges used. For this reason the algorithm of [BJ, to the best of our knowledge, can not determine 
the parity of the number of all Hamiltonian cycles. 

Recall that the algorithm of [H O [23] uses dynamic programming based on the observation 
that if we have constructed a path from v a to Vb, then the only additional information needed 
is the set of vertices used along the way. Hence the essential information required to store for 
all subpaths is their endpoints and the set of visited vertices, leading to 0(n2 n ) table entries. 
Recall that for designing a fast dynamic programming algorithm, it is essential that one chooses 
the right strategy of decomposing a candidate solution. For example, if G has a small balanced 
separatoiQ of size 0.1 |V| it seems impossible to modify the above algorithm in order to exploit 
this fact. 

In this work we study the following way to decompose a Hamiltonian cycle H: Consider 
an arbitrary ordering ei,...,e m of the edges E and decompose H for every i into H\ = 
{e%, . . . , ei} D H and H2 = {ej+i, . . . , e n } n H. This suggests another dynamic programming 
strategy alternative to [H El [23] that was much less well-studied, mainly because it seemed a 
priori that the essential information needed to store for a set of partial solutions H\ is 2 e ( nlgn ): 
The degrees (zero, one, or two) of each vertex with respect to H\ and the pairing in which 
vertices of degree one are connected by H\. The informed reader will notice that this way of 
decomposing solutions is inherent to dynamic programming algorithms that operate on tree 
decompositions or path decompositions. 

The notions of pathwidth / treewidth proved to be an excellent tool for dealing with many NP- 
hard problems on graphs. In the 1970s and 1980s, several groups of researchers discovered the 
concept independently. In their fundamental work on graph minors, Robertson and Seymour 
|37| introduced the notions path /treewidth and path/tree decomposition, and these became the 
dominant terminology. A graph having small path/treewidth means informally that it can be 
decomposed efficiently in a path/tree- like manner. Many problems can be solved using dynamic 
programming by decomposing a solution according to the given path/tree decomposition. We 
refer to [8] for more information on these notions. 

In early work, such as the influential result of Courcelle [13], algorithms with running times 
of the type /(pw)n were giverj^J Later it was noticed that for many such algorithms the function 
/(pw) can be substantially improved, greatly improving the tractability as well. For example for 
many problems whose solutions can be verified by separately considering its intersection with all 

1 A vertex set X such after removing X, all connected components have size at most |V|/2. 
2 Assuming an input graph G on n vertices along with a path decomposition of G width pw is given. 
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neighborhoods of the vertices of the input graph (see [22] ) , we can obtain 2°^n oi ^> by employ- 
ing dynamic programming. For example the Independent Set, or equivalently, Vertex Cover, 
problem can be solved in 2 pw n time [28, 34] . From the work of [25J and standard Karp-reduction 
it follows that this dependence cannot be improved to sub-exponential algorithms unless the 
Exponential Time Hypothesis fails, i.e., unless CNF-SAT has a subexponential algorithm. In 
|32| it was shown that under a stronger assumption (the so-called Strong Exponential Time Hy- 
pothesis (SETH)), the current algorithms are optimal even with respect to polynomial jumps, 
that is, problems with current best running time f{tw)n 0<y1 ^ cannot be solved in /(tw) 1 ~ <E n°( 1 ) 
for positive e where /(tw) is 2 tw , 3 tw for respectively Independent Set and Dominating Set. 

For the Traveling Salesman problem, or even the Hamiltonian cycle problem, the situa- 
tion is different since they are not local problems (see for example |221 Section 5]). In |16j . 
the notion of a Catalan structure was introduced to give an algorithm that solves the /c-path 
problems on f/-minor free graphs in 2 ( tw ) time. To obtain this, their result works for 'pairing- 
encodable' problems (i.e. problems where the connectivity properties can be encoded by a 
matching/pairing) since it bounds the number of ways paths can intersect with a part of the 
tree decomposition. However, in [15] it was shown that along with many other connectivity 
problems, Hamiltonian cycle can be solved by a Monte Carlo algorithm in 4 c> ( tw )n C) ( 1 * ) time. 
For many of these algorithms with running time /(tw)n c,( - 1 - ) , it was shown that an algorithm 
running in time /(tw) 1_<E n c '( 1 ) would violate SETH. Nevertheless, the exact complexity of the 
Hamiltonian cycle problem remained elusive. 

Recently, a superset of the current authors, found a connection of the optimal substructure 
of a dynamic programming algorithm with the rank of a certain matrix [9]: Consider a matrix 
H with rows and columns indexed by partial solutions, with a 1 if and only if the two partial 
solutions combined give a valid solution, then if we have more than rk(%) partial solutions, one 
will be redundant in the sense that we can safely forget it. This leads to deterministic 2°^ tv ^n c> ^ 
time algorithms for many connectivity problems, in a sense derandomizing the work of Cygan 
et al. |15j . but the used techniques do not match their runtimes, e.g., the ^'n ' 1 ' time for 
Hamiltonian cycle. 



Our contribution 

Inspired by the mentioned result from [9], we study a matrix that we call the Matching Con- 
nectivity matrix. We present a family of perfect matchings Xj, which we show to be a basis 
of the Matching Connectivity matrix, therefore we establish its rank. All further results of the 
paper use the basis Xj and its properties as its key tool, which we now elaborate on. 



Determining rank of the Matching Connectivity matrix. For an even integer t > 2, 
the Matching Connectivity matrix Tit is a matrix that has rows and columns both labeled 
by all perfect matchings of the complete graph Kt on t vertices; an entry Ht[Mi, M2] is 1 
if Mi U M2 is a Hamiltonian cycle and otherwise. The centerpiece of our work is that Tit 
has rank exactly 2*/ 2-1 over GF(2). To establish this result, we define an explicit family X^ 
of 2*/ 2-1 perfect matchings on Kt and show that its columns (or rows) form a basis of Tit over 
GF(2): First, essentially by design, each matching M £ Xj has a unique partner M' E Xj such 
that their union is a Hamiltonian cycle. Thus, the rows and columns labeled by Xf induce a 
permutation matrix, which implies the required rank lower bound. Second, we give an explicit 
factorization of Tit into a product of two rectangular matrices with inner dimensions indexed 



by ~Kt in Theorem 3.4 This proves that Xf is indeed a basis, provides an explicit formula for 
linear combinations, and completes the claimed rank bound (see Section [3]). 

In |36j a matrix Ti' t is studied that is obtained by restricting Tit to all perfect matchings of the 
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complete bipartite graph with independent sets of size t/2. There it is shown that when taken 
over the reals, the rank of %[ is exactly ( t *2-i) using group representation theory. Using this 
the authors disprove the originaj^] 'log-rank conjecture' in communication complexity since they 
also showed that the non-deterministic communication complexity of determining whether two 
such matchings give a Hamiltonian cycle is O(nloglogn). In [9j Lemma 3.13] a factorization of 
a matrix that contains Tit as submatrix into two matrices with inner dimension 2* _1 was given. 



Exact algorithms for Hamiltonicity. By exploiting the peculiar form of the basis X n we 
show that the number of distinct subsets of all the matchings of X ra is C(1.888 n ). Together 
with Theorem 3.4 this allows us to show deterministic algorithms computing the parity of the 
number of Hamiltonian cycles in undirected graphs and directed bipartite graphs (section 3.2.2) 
in 1.888 n rt°( 1 '' time. By combining those results with the Isolation Lemma we obtain Monte 
Carlo algorithms solving the decision version of Hamiltonicity in both undirected and directed 
bipartite graphs within the same running time. 

Even though our algorithm for undirected graphs is slower than the algorithm of Bjorklund |6j , 
we believe it is of interest as it uses very much different tools and allows solving the problem 
also in directed bipartite graphs. We would like to recall that solving the Hamiltonicity problem 
on undirected bipartite graphs in 0{{2 — e) n ) time was the first step of Bjorklund on the way to 
the algorithm for general undirected graphs. For this reason we believe that studying directed 
bipartite graphs is justified. 



Algorithm for bounded pathwidth. Using the set of perfect matchings X( we obtain a 
faster algorithm for solving the Hamiltonian cycle problem via a non-trivial pathwidth dynamic 
programming routine. Regarding the algorithm, already from the rank bound it can be argued 
that only (2 + \/2) pw space is needed to solve Hamiltonian cycle on a graph with given path 
decomposition of width pw. The key idea is to replace memoization of all partial solutions 
by storing only fingerprints of groups of solutions that encode how many partial solutions are 
consistent with a given basis matching. In fact, these numbers are stored only modulo two 
which requires only one bit per basis matching. To achieve also time (2 + \/2) pw n c '( 1 ^ we show 
in Lemma 4.1 how to efficiently convert these fingerprints from one basis to another, which 
permits us to perform the needed dynamic programming table computations (in particular 
insertion of edges into partial solutions); this crucially depends on the structure of X 4 . Notably, 
Lemma 4.1 gives a second proof of the rank upper bound for Tit, but it does so in a more 
implicit way; in particular it does not provide an explicit factorization (see Section [4]). 

Let us point out the main differences to the related algorithmic results obtained by Bod- 
laender et al. [9]: The present faster algorithm for Hamiltonicity parameterized by pathwidth 
required a new dynamic programming strategy, unlike the results of [§] that speed up existing 
formulations. Furthermore, we require randomization (for the Isolation Lemma |33j) to guar- 
antee a unique solution, in order for the fingerprinting approach to work. Finally, achieving 
time (2 + \/2) w n c ' ( - 1 ^ depends crucially on the structure of our basis matchings X^ (to allow for 
Lemma 4.1) and does not appear to follow directly from the rank. (Similarly, the subsequent 
lower bound requires the existence of a sufficiently large permutation/identity matrix in Tit and 
does not follow directly from the rank lower bound.) 



Matching lower bound assuming SETH. We show that if the running time of our algo- 
rithm can be significantly improved, then satisfiability of CNF-Sat formula's of m clauses and n 

3 A weakened version is still a standing open problem in communication complexity [2l Section 13.2]. 
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vertices can be determined in (2 — e) n mP^ time. The latter would contradict the Strong Expo- 
nential Time Hypothesis (SETH) introduced in [25]. Although there is not consensus about it's 
truth, a number of results have been given assuming SETH [35, IU[32j[T5]. As with previous 
results, this result should be interpreted that there is a boundary to significantly improving our 
algorithm, namely finding the so far elusive (2 — e) n mP^ algorithm for CNF-Sat. 

For the tight runtime lower bound we use our basis as a part of a gadget in a reduction from 
CNF-Sat. Although the basic setup is similar to previous lower bounds [32j [15] our reduction 
is different in the sense that we require a very generic gadget (the induced subgraph gadget, 
discussed in Subsection 5.1 ). Using this, we can exploit the crucial property that the submatrix 
of the Matching Connectivity matrix that is induced by the columns and rows of Xt is a 
permutation matrix, i.e., each basis matching has a unique partner in Xt such that their union 
is a Hamiltonian cycle. Then, in a similar but technically challenging vein to the reductions 
of [32l [15] , choices can be transferred through a series of gadgets (see Section [5]) . 



Further algorith mic c onclusions. As a corollary of our bounded pathwidth algorithm in 
Section |Z| (Corollary [iTE]) we also obtain a (2 + v / 2) ?1 / 6 n c '( 1 ) = 1.1583 n n°W time Monte Carlo 
algorithm for Hamiltonicity in cubic graphs, which to the best of our knowledge is the fastest 
known algorithm in this class of graphs. 



2 Preliminaries 

Graphs. We use standard graph notation. For a graph G = (V, E) we denote V(G) and E{G) 
for its vertex and edge set respectively. For X, Y C V we let E(X, Y) be the set of all edges 
with one endpoint in X and one in F. A Hamiltonian cycle is the edge set of a simple cycle 
that visits each vertex exactly once. A cycle cover is a set of edges F C E such that each vertex 
of G is incident with at exactly two of these edges (i.e., the edges form cycles). In a partial cycle 
cover each vertex is incident with at most two edges (i.e., the edges form paths and cycles). 



Perfect matchings. A perfect matching of a graph is a set of edges such that each vertex is 
incident with exactly one of them. It is well known that the union of any two perfect matchings 
in a graph forms a cycle cover of the graph (where some cycles are potentially of length 2). 
Given some base set U, we use Il2(i7) for the set of all perfect matchings of U (i.e., if U has no 
graph structure then all partitions into sets of size two each is included). Borrowing from the 
partition lattice partially ordered by refinement, we use M\ n Mi = {U}, for M±, Mi € U2(U), 
to express the fact that the union of the two perfect matchings Mi and M2 is a Hamiltonian 
cycle; for two perfect matchings this is equivalent to getting the trivial partition {U} into a 
single set as the outcome of the meet-operation n . We do not require any further tools or 
notation from the partition lattice. 



Pathwidth and path decompositions. A path decomposition of a graph G = (V, E) is 
a path P in which each node x has an associated set of vertices B x C V (called a bag) such 
that (J B x = V and the following properties hold: 

1. For each edge {u, v} £ E{G) there is a node x in P such that u, v £ B x . 

2. If v & B x n By then v £ B z for all nodes z on the (unique) path from x to y in P. 

The pathwidth of P is the size of the largest bag minus one, and the pathwidth of a graph G 
is the minimum pathwidth over all possible path decompositions of G. Since our focus here 
is on dynamic programming over a path decomposition we only mention in passing that the 
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related notion of treewidth can be denned in the same way, except for letting the nodes of the 
decomposition form a tree instead of a path. 

It is common for the presentation of dynamic programming to use path and tree decom- 
positions that adhere to some simplifying properties, in order to make the description easier 
to follow. The most commonly used notion is that of a nice tree decomposition, introduced 
by Kloks [29J; the main idea is that adjacent nodes can be assumed to have bags differing by 
at most one vertex (this can be achieved without increasing the treewidth). For an overview 
of tree decompositions and dynamic programming on tree decompositions see |1CH I24j . In a 
similar way, but using also the extension of introduce edge bags from |15] we define nice path 
decompositions as follows. 

Definition 2.1 (Nice Path Decomposition). A nice path decomposition is a path decomposition 
where the underlying path of nodes is ordered from left to right (the predecessor of any node is 
its left neighbor) and in which each bag is of one of the following types: 

• First (leftmost) bag: the bag associated with the leftmost node x is empty, B x = 0. 

• Introduce vertex bag: an internal node x of P with predecessor y such that B x = 
By U {v} for some v ^ B y . This bag is said to introduce v. 

• Introduce edge bag: an internal node x of P labeled with an edge {u, v} G E{G) with 
one predecessor y for which u,v G B x = B y . This bag is said to introduce uv. 

• Forget bag: an internal node x of P with one predecessor y for which B x = B y \ {v} for 
some v £ By. This bag is said to forget v. 

• Last (rightmost) bag: the bag associated with the rightmost node x is empty, B x = 0. 

It is easy to verify that any given path decomposition of pathwidth pw can be transformed 
in time iV^G^Ipw^ 1 ) into a nice path decomposition without increasing the width. 



Further notation. For two integers a, b we use a = b to indicate that a is even if and only if 
b is even. We use Iverson's bracket notation: if p is a predicate we let [p] be 1 if p if true and 
otherwise. If u : XJ — > {1, . . . ,iV}, we shorthand uj(S) = Ylie&S u i e ) f° r S C U. 



3 Structure of the Matching Connectivity matrix and applica- 
tions 



The section is outlined as follows: We will first determine the structure of the Matching Con- 



nectivity matrix Tit in Subsection 3.1 More specifically, we give a basis and determine the rank 



exactly through an explicit matrix factorization. In Subsection 3.2 we will give an application 
of the matrix factorization to exact exponential algorithms for Hamiltonicity. 



3.1 Bases and factorization of Matching Connectivity matrix. 

Recall that the matrix Tit has rows and columns both labeled by all perfect matchings of the 
complete graph K-t on t vertices; an entry Ht[Mi, M2} is 1 if M\ U M2 is a Hamiltonian cycle 
and otherwise. The dimension of the matrix is 

t\ t\ 
(i)!-2f X (i)!-2l' 
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Figure 1: The matrix T~Lq. Letting the baseset be {0, . . . , 5} matching 1 indexing row and column 
1 equals {{0, 1}, {2, 3}, {4, 5}}. The set X 4 = {1, 2, 4, 5} from Definition 3.1 is easily seen to be 
a row basis; the linear combinations are depicted in the last column. 



Let us point our some small cases: For t = 2 we have only one perfect matching on the two 
vertices, and with the union of two such matchings being considered as a Hamiltonian cycle (in 
all other cases, where t > 4, there cannot be a Hamiltonian cycle if the two perfect matchings 
have at least one edge in common). For t = 4, there are 3 perfect matchings and H3 is easily 
seen to be the complement of the 3x3 identity matrix. The matrix is a 15 x 15 matrix 
shown in Figure [TJ 

To prove the exact value of the rank we introduce for each even t > 2 a family Xt of 2*/ 2-1 
perfect matchings with the goal of proving that the corresponding columns (or rows) of T~Lt form 
a basis over GF(2). The definition of Xj requires the vertices to be ordered, say 0, 1, . . . , t — 1, 
and edges in the matchings are very local, i.e., their endpoints are at most at distance three 
with respect to the ordering. From the structure of the matchings it will be easy to see that 
they give a lower bound for the rank of %t ■ The submatrix of %t induced by rows and columns 
from X( is a permutation matrix, which already has rank 2*/ 2-1 itself. This property will be of 
the essence for our lower bound on the runtime of pathwidth-based dynamic programming for 
the Hamiltonian cycle problem in Section [5] 

Getting the matching upper bound is more involved. We obtain this result by giving a 
concrete factorization of Tit in terms of two rectangular submatrices of %t induced by the rows 
respectively columns X< in %t] this is done by a rather technical inductive argumentation. This 



will form the basis of our algorithm from Subsection 3.2 



Let us begin by introducing the families X^. For a perfect matching M S ^(i/) we define 
a function om- U — )• U with «a/W = 3 if and only if {i,j} E M, i.e., maps each element 
of U to its partner in the perfect matching M. We shorthand Ut ■= {0, 1, . . . , t — 1}. 
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Definition 3.1. Let e denote the empty string. We let X(2, e) := {{0, 1}} and X2 := {^(2, e)}. 
Let t > 4 be an even integer and let a be a bit-string of length | — 2. We define perfect 
matchings X(t, a0) and X(i, al) of ^ = {0, . . . , t — 1} as follows: 

X(t, al) := X(t - 2, a) U {{t - 2, t - 1}}, 

X(i, a0) := (X(t - 2, a) \ {{* - 3, a(/ - 3)}}) U {{i - 2, a(t - 3)}, {t-3,t- 1}}. 

We shorthand X(a) for X(2|a| + 2, a) since the bitstring a determines the size t of the base set 
(namely Ut). We use a to denote the binary complement of a bit-string a. Finally, we let Xt 
be the set of all perfect matchings X := X(t,a) for any bitstring a of length | — 1. 

Unfortunately, the formal definition is not very enlightening regarding the actual structure of 
the perfect matchings X(t, a) in X<. Let us clarify their structure by first showing the matchings 
for t = 4 and t = 6; recall that X(2, e) = {{0, 1}}. 

X(4, 1) = {{0, 1}, {2, 3}} X(4, 0) = {{0, 2}, {1, 3}} 

X(6, 11) = {{0, 1}, {2, 3}, {4, 5}} X(6, 10) = {{0, 1}, {2, 4}, {3, 5}} 

X(6, 01) = {{0,2}, {1,3}, {4, 5}} X(6,00) = {{0, 2}, {1, 4}, {3, 5}} 

When recursively constructing further perfect matchings for some even integer t > 4, we always 
either add another edge on the two new elements (when the new bit is 1) or we replace the last 
edge (i.e., the one matching t — 3 to t — 4 or t — 5; note that i — 3 is the last element for t' = t — 2) 
by matching its vertices to the two new elements (when the new bit is 0). 

Let us explain the intuition behind the bitstrings: Let t = 6 and group the elements as | 
1,2 I 3,4 I 5. Observe that the matchings X(t, ■) are exactly all choices of matching the elements 
such that each edge connects two elements that have exactly one dividing vertical line between 
them, i.e., all choices of perfect matchings that match only elements from adjacent groups. Now 
the first bit in the bitstring determines whether the first edge is {0, 1} or {0, 2} (these are all 
possible options for matching under the group restriction). Depending on this either 2 or 1 
still needs to be matched to 3 or 4; the latter choice is determined by the second bit. The last 
edge must always go to element t — 1 (i.e., 5 in this example), so there are only | — 1 bits. 

Proposition 3.2. Let t > 2 be an even integer, and group the elements {0, ... ,t — 1} into | 
1, 2 I . . . I i — 3, t — 2 \ t — 1. The family X t consists of all perfect matchings that match only 
elements from different, but adjacent groups. There are 2 t l 2 ~ 1 such perfect matchings. 

Clearly, the presented families of perfect matchings, one for each even integer t, have a very 
particular and symmetric structure. Our aim is to show that the 2*/ 2 " 1 perfect matchings form 
a basis for the Matching Connectivity matrix %f Regarding any two such matchings X(t,a) 
and X(t, b) it is not hard to show that their union is a Hamiltonian cycle if and only if a = b. 

Proposition 3.3. Let t > 4 be an even integer and let a and b be bit-strings of length | — 1. 
Then X(t,a) U X(t,b) is a Hamiltonian cycle of K t , or equivalently, X(t,a) n X(t,b) = {Ut}, 
if and only if b = a. 

Proof. Consider the first position, say i, such that a[i] = b[i] (we consider a and b to be indexed 
from left to right, starting with 0). All earlier positions j < i are hence different, and following 
the definition of X(t, ■) this means that they prescribe exactly opposite choices. E.g., one 
matching will match to 1 and the other matches it to 2, without loss of generality {0, 1} G X(a) 
and {0,2} G X(b). Consequently, the next bit specifies the matching for 2 in X(a) and the 
matching for 1 in X(b). This pattern continues and effectively we obtain two alternating paths 
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that start from and follow edges from X(a) and X(b) alternatingly. If the paths meet only 
at t — 1 (for which there is no bit that allows an alternate choice of matching) then together 
they give a Hamiltonian cycle. Since we assumed that a[i] = b[i], the paths meet when the 
bits prescribe that both matchings match to the same element (it can be verified that bit i 
decides whether the so far unmatched element of 2i — 1 and 2i is matched to 2i + 1 or 2i + 2). 
Thus, we have found a closed cycle in the union of X(a) and X(b) which does not contain all 
vertices, hence they do not form a Hamiltonian cycle. (Note the special case of a[0] = b[0] which 
indicates that both matchings contain {0,p} for p £ {1, 2}.) □ 



From Proposition 3.3 we directly get a rank lower bound of 2 l l 2 1 since the submatrix given 
by all rows and columns of matchings X(t,-) is a permutation matrix of size 

2 t/2-l x 2 */2-l. 

Now we state the main theorem of this section. Due to its technicality, the proof is deferred to 
Section [A} The statement of the theorem is equivalent to saying that the Matching Connectivity 
matrix Tit can be written as the product of two rectangular submatrices of %t whose rows 
respectively columns are labeled by matchings from X^. This implies that the set of those 
rows/columns forms a basis for Ht and that its rank is 2*/ 2-1 . 

Theorem 3.4. Let t>2 be an even integer and let Mx,M% £ ~n.2(Ut). It holds that 

[Mi n M 2 = {U t }] = Yl t Ml n x ^ a ) = Ml • [ M2 n x (*>") = 
ae{Q,iy/ 2 - 1 



where X(t,a),X(t,a) E Xj according to Definition 3.1, (Each matching in Xj occurs exactly 
twice, once as X(t,a) and once as X(t,a).) 

Corollary 3.5. The rank of the Matching Connectivity matrix Tit over GF(2) is 2'/ 2 - 1 for all 
even integer t>2. 



Proof. It follows from Proposition 3.3 that the rank is at least 2*/ 2-1 : The Matching Connec- 
tivity matrix Tit contains a 2*/ 2 ~ 1 x 2*/ 2-1 submatrix induced by the columns and rows of all 
perfect matchings in Xf which is a permutation matrix. 



From Theorem 3.4 we immediately get that the rank is at most 2*/ 2_1 : We can read the 
theorem statement as a factorization of Ht as the product of two submatrices of %f The first 
submatrix has rows labeled by all perfect matchings on Ut and columns labeled by basis match- 
ings X(t,a) for lexicographically ordered bitstrings a. The second matrix has columns labeled 
by all perfect matchings of Ut and rows labeled by basis matchings X(t,a) for lexicographi- 
cally ordered bitstrings a. The rank upper bound follows immediately from the fact that both 
matchings have rank at most 2*/ 2-1 (corresponding to their smaller dimension). □ 

3.2 Exact exponential algorithms for Hamiltonicity 

In this section we present Monte Carlo algorithms for solving the Hamiltonian cycle problem in 
time 0(1.888 n poly(n)) in undirected graphs and directed bipartite graphs. These algorithms 
are based on further ideas and insights about the families Xt of perfect matchings, and in 



particular we greatly rely on Theorem 3.4 



First, we show that to solve the decision version it is enough to solve the problem of computing 
the parity of the number of Hamiltonian cycles modulo two. The main part of our algorithm 



lies in the proofs of the following two lemmas (the proofs are provided in Section 3.2.2) 



Lemma 3.6. There is an algorithm, which given an undirected graph G = (V, E) together with 
a weight function oj : E — > {1, . . . ,u> max } finds the parity of the number of Hamiltonian cycles 
of weight w for every w G [0, . . . , nw max ] in C(1.888 n poly(n + w max )) time. 
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Lemma 3.7. There is an algorithm, which given a directed bipartite graph G = (V,A) together 
with a weight function uj : A — > {1, ... ,aj max } finds the parity of the number of Hamiltonian 
cycles of weight w for every w G [0, . . . , nu max ] in 0(1.888 n poly(ra + u> max )) time. 

Now, by an application of the Isolation Lemma, we can show that our modulo two counting 
of solutions suffices to determine (with high probability) whether or not G is Hamiltonian. 

Definition 3.8. A function lo : U — > Z isolates a set family T C 2 U if there is a unique S' G T 
with uj(S') = minsgjr oj(S). 

Lemma 3.9 (Isolation Lemma, [33J). Let J- C 2 U be a set family over a universe U with 
\J-\ > 0. For each u £ U , choose a weight u)(u) £ {1, 2, . . . , N} uniformly and independently at 
random. Then probfcj isolates J-] > 1 — \U\/N . 

Theorem 3.10. There exists a Monte Carlo algorithm solving the Hamiltonian cycle problem 
in 0(1.888 n poly(ra)) time in undirected graphs and directed bipartite graphs. 

Proof. Given a graph G (either undirected, or directed bipartite) with m edges (arcs), for each 
edge (arc) assign an integer weight from the interval [1, . . . ,2m] uniformly and independently 



at random. Then we use Lemma 3.6 (Lemma |3.7| ) to calculate the parity of the number of 
Hamiltonian cycles of each weight w £ [0, . . . , nw max ]. If for some w there is an odd number of 
Hamiltonian cycles, then our algorithm returns YES, otherwise it returns NO. 

The running time of the algorithm follows from the running time of the black-box usage 
of the parity calculating algorithm. If there is no Hamiltonian cycle in our graph, then our 
algorithm certainly returns NO. However, if the graph contains at least one Hamiltonian cy- 



cle, then by Lemma 3.9 with probability at least 1/2 our weight function isolates the family 
of all Hamiltonian cycles of G and consequently for some weight there is an odd number of 
Hamiltonian cycles and our algorithm returns YES. Therefore we have obtained a Monte Carlo 
algorithm. □ 

3.2.1 Further uses of the basis matchings 

In this section we give two technical lemmas that form the core of our two algorithms, which 
are based on the families X n of perfect matchings introduced in Section [3} First, we show 



that the number of subsets of all matchings in X n is bounded by C(1.888 n ) (Lemma 3.11) 



Second, we show how to compute the number of extensions of basis matchings to Hamiltonian 



cycles (Lemma 3.14); for this we use dynamic programming over the mentioned subsets of basis 
matchings. 

Lemma 3.11. The set of all subsets of all matchings from the basis X n contains C(1.888 n ) 
distinct matchings. 

Proof. Recall that X n = {X(n,a) \ a 6 {0, l}™/ 2 " 1 }. For an even integer n define: 
t( n ) = \{S C X(n, o) | o € {0, l}"/ 2 - 1 }] 

t 2 { n ) = \{S C X(n + 2, a) | a G {0, l} n/2 A S n {{n, n + 1}, {n - 1, n + 1}} = 0}| . 

Less formally, t(n) is the number of distinct matchings being subsets of the basis, whereas for 
t^ip) we consider longer bitstrings, and additionally we assume that the subset does not contain 
an edge incident to the last element of the universe, i.e., to n + 1. In what follows we prove the 
following inequalities: 

t(n) < 2t(n - 2)+t 2 {n - 2) , 
t 2 (n) < 4i(n - 2) + t 2 (n - 2) . 
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The first inequality follows from the case analysis of the last bit of the string a and whether 
S contains the edge incident with n — 1 or not (which we refer to as the last edge). When we 
analyze £(n) we have 4 cases: 

• the last bit of a is 0, the set S contains the last edge, 

• the last bit of a is 1, the set S contains the last edge, 

• the last bit of a is 0, the set S does not contain the last edge, 

• the last bit of a is 1, the set S does not contain the last edge. 

Note that all the subsets S from the first two cases can be upper bounded by 2t(n — 2), and the 
from the second two cases by t2(n — 2), which follows directly from the definition of ti- Now we 
analyze t2(n), here we have 8 cases to consider, depending on the last two bits of a (4 choices), 
and whether S contains the edge between {n — 3, n — 2} and {n — l,n} (2 choices), which we 
call the penultimate edge (note that the edge incident to n + 1 is definitely not contained in S) . 
When S contains the penultimate edge, then we bound each of the 4 cases depending on the 
last two bits of a independently by t(n — 2). Observe that so far we made no savings and no 
cases we considered identical. However the key case is when S does not contain the penultimate 
edge, which means that S contains neither the penultimate nor the last edge. All 4 such cases 
can be upper bounded by ti(n — 2), as the last bit of a does not matter anymore. 

Observe that £(2) = 2 and £2(2) = 3. Moreover if we multiply the horizontal vector 

(£(ra), £2(71)) by the matrix A = I j we obtain (a, b) where a > t(n + 2) and b > £2(71 + 2). 

At the same time A = BDB -1 , where 




1 




Vu 


2VT7 


1 






2\fY7 




Consequently A n = BD n B~ 1 where 



3+ 2 ) ) an d the claimed upper bound follows, since matrices B and 
B^ 1 contain only fixed constants which are hidden inside the O-notation. In fact one can show 
that our analysis is tight, since we do not overcount any subsets S neither in £(n) nor in £2(71) 
in our case analysis. □ 

Definition 3.12 (perfect matchings). For an undirected graph H by I^-H") we detote the 
set of perfect matchings in H. 

Recall that by 112(F) we denote the set of all matchings in the complete graph on V. 

Definition 3.13 (extensions to Hamiltonian cycle). For an undirected graph H = (V,E), 
a weight function uj : E — > {0, . . . , w max }, a matching M G Il2(V) and an integer w by 
ext(M, w, H, 00) we denote the number of perfect matchings in H of weight w which together 
with M form a Hamiltonian cycle, i.e., 

ext(M,w,H,u)= Yl [M'nM = {V}]. 

M'eU 2 (H),u>(M')=w 
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Lemma 3.14. Given an undirected graph G = (V, E) and a weight function uo : E — )• {0, . . . , w max } 
one can in 0(1.888 n poly(n+w max )) time compute all the (na; ma x+l)2 n / 2_1 values ext(X(n, a),w, G, uj) 
for a G {0, l} n / 2_1 and < w < nw max . 

Proof. For a (not necessarily perfect) matching M in G, a vertex u and weight < uu < noj max 
let t[M] [v] [w] be the number of walks in the complete graph on V from an arbitrary fixed vertex 
v\ to the vertex v which (i) contain an even number of edges, (ii) start with an edge of M, (iii) 
alternately contain an edge of M and an edge of E, (iv) use each edge of M exactly once (v) the 
total weight of edges from E equals w. Observe that ext(M,w,G,ui) = t[M][v{\[w], hence it is 
enough to describe a dynamic programming routine computing all the values t[X(n, a)] [vi] [w] 
for a G {0, l}™/ 2 " 1 , < w < nw max . 
We use the recursive formula 

t[M \ {uu}} [u] [w - u(u, v)] . 

uu'&M 
u£N(v) 



We also define a corner case t[0][i>][iu] = [v = v i A w = 0] for v G V. By Lemma 3.11 the above 



formulas together with memoization proved the claimed algorithm. □ 
3.2.2 Proofs of Lemmas [376] and [3771 



In this section we prove Lemmas 3.6 and 3.7, that is we focus on computing the number of 



weighted Hamiltonian cycles modulo two. 



Proof of Lemma 3.6, First, we iterate over pairs of distinct edges e±, e2 incident to some arbi- 
trary fixed vertex v\. We want to find the parity of the number of Hamiltonian cycles containing 
both e\ and e%. Note that for a different pair of edges incident to v\ we count different Hamil- 
tonian cycles and each Hamiltonian cycle is counted exactly once. Since we need an additional 
assumption that n is even, if this is not the case we subdivide the edge e%, resulting in two new 
ed ges e'j, and we set cj(e^) — w(ei), ui{e!^) — 0, where is incident to v%. We somewhat 
abuse notation and by e\ denote the edge e[. Observe that there is a bijection between weighted 
Hamiltonian cycles before the subdivision and after the subdivision. 

Note that any Hamiltonian cycle containing e\ and e2 can be uniquely decomposed into two 
perfect matchings Mi, Mi in G where e\ G M\ and ei G Mi- Let us fix a weight < w < nui maJi . 
Our goal is to compute the number of pairs {M\,M2), where M\,Mi are perfect matchings in 
G, ei G Mi, e% G M2, u(Mi) + u(M%) = w and Mi U M2 is a Hamiltonian cycle. Formally, we 
want to calculate the following sum modulo two 

£ [Mi n M 2 = {V}} . 

Mien 2 (G) M 2 en 2 (G) 

eiSMi e 2 eM 2 

u(M-l)+u(M 2 )=w 
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By Theorem |3.4| the above formula is equivalent (modulo two) to 

E E E [M 1 nX(n,a) = {V})-[M 2 nX(n,a) = {V}] 



Mit=n 2 (G) M 2 en 2 (G) aeio.i}"/ 2 - 1 

eiEMi e 2 SM 2 

uj{M{]+u{M 2 )=w 

= E E [MinX(n,a) = {y}] £ [M 2 nX(n,a) = {y}] 

ae-to.i}"/ 2 - 1 A/ien 2 (G) M 2 en 2 (G) 

eieMi e 2 eM 2 

= y~] ext(X(ra, a),u>i, (V, • ext(X(n,a), w 2 , (V", E 2 ),u). 

aglO,!}™/ 2 " 1 wi+w 2 =i« 

In the above, we first change the summation order, and then for i = 1,2 denote Ei = (E \ 
E(vi,V)) U {e^}. Note that all perfect matchings contributing to ext(X(n, a), E\), to) 
(and analogously for E 2 ) contain the edge ei, as this is the only edge in the graph (V, E\) 
incident to v\. 

By Lemma 3.14 we can find all the needed values in the claimed running time, hence the 
proof of Lemma 3.6 follows. □ 

In the proof of Lemma 3.7 we have to deal with directed graphs, while still using Lemma 3.14| 
which can only handle undirected graphs. However as we will show one can exploit the bipar- 
titeness to harness the directedness of the graph. 



Proof of Lemma 3. 7. Let G = (V± l±) V 2 , A) be a directed bipartite graph. Clearly we can assume 
l^i I = | V 2 \, since otherwise there are no Hamiltonian cycles in G. We create two auxiliary 
weighted undirected bipartite graphs Gt = (V\ l±) V 2 ,Ee,uJi) and G r = (V\ ttJ V 2 , E r ,oj r ) where 
Ei = {uv : u G Vi,v 6 V 2 , (v,u) £ A}, E r = {uv : u G Vi,v G V 2 ,(u,v) £ A} and for u £ V±, 
v G V 2 we have uj r (uv) = uj(u, v) and loi(uv ) = uj(v, u). That is we split the arcs of A depending 
on whether they have their start-point in V\ or V 2 and take the two underlying undirected 
graphs after the split. 

Note that each Hamiltonian cycle in G can be uniquely split into two sets of arcs, one of 
which corresponds to a perfect matching in Gi and the other in G r . Moreover any pair of a 
perfect matching in Gi and a perfect matching in G r together forms a cycle cover in G, which 
might consists of several cycles. Our goal is to consider all pairs of perfect matchings in Gi and 
G r which together form a Hamiltonian cycle in the underlying undirected graph of G, which 
guarantees that we count exactly the Hamiltonian cycles in the directed bipartite graph G. Let 
us fix an integer < w < nw max and count the number of Hamiltonian cycles of weight w in 
the graph G. 

E [Mi n M 2 = {V}] 

Mien 2 (Gf) M 2 en 2 (G r ) 

Lu e (Mi)Uuj r (M 2 )=w 

(Thm\M) = E E E t Ml n X ( n > a ) = " [ M2 n X ( n ^) = 

Af!en 2 (Gf) M 2 en 2 (G r ) ae{o,i}™/2-i 

LJl(Ml)UL0 r (M2)=W 

= E E [MinX(n,a) = {V}) £ [M 2 n X(n,a) = {V}] 

ae{o,i}«/ 2 - 1 M 1 en 2 (G^) M 2 en 2 (G r ) 

aj ( (Mi)Uw r (M 2 )=w 

= ext(X(n,a),wi,Gi,u;£) ■ exi(X(n,a),w 2 ,G r ,uj r ). 

aefO,!}™/ 2 - 1 w 1 +w 2 =w 



By Lemma 3.14 the proof of Lemma 3.7 follows. □ 
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4 Solving Hamiltonian cycle fast on pathdecompositions 



In this section we present an re 3 (2 + \/2) pw pw c, ^ :L ^-time algorithm for solving Hamiltonian cycle 
on a graph G with a given path decomposition of width pw. Recall that partial solutions for 
Hamiltonian cycle are sets of paths such that all vertices before the current bag are internal 
in some path, and vertices in the current bag may be endpoints, internal, or unused. It then 
suffices to remember for each such partition of the current bag, in what way the endpoints are 
connected into pairs (these arrangements are perfect matchings on the set of endpoints); it is 
well-known that any further information about the paths is not needed. The downside is that 
this involves roughly pw pw many partial solutions which dominates the runtime. 

The key idea for our much faster algorithm is as follows: Instead of storing for all partitions 
into endpoints, internal, and unused vertices all the possible perfect matchings of the endpoints, 
we only store, intuitively, a combined "fingerprint" of all matchings together. Indeed, we fix 
an ordering of the vertices and store for each matching of the resulting family Xt the number 
of partial solutions that give a single cycle together with this matching. (These matchings 
abstract away the need of connecting through all so far unused vertices since this is covered 
by the partitions.) In fact, since our basis works only over GF(2) we count those solutions 
modulo two. This however still useful since we can essentially ensure the existence of a unique 
Hamiltonian cycle of minimum weight via the Isolation Lemma (and we need to solve a weighted 
version of our modulo two counting problem) . 

Given this setup, let us solve the following problem by dynamic programming on a path 
decomposition: Given a graph G = (V, E) along with a path decomposition of pathwidth pw, 
and non- negative edge weights u: E — > {1, . . . ,w max }- The task is to compute for each ui* G 
{1, . . . ,n ■ uj max } the parity of the number of Hamiltonian cycles of G with weight exactly u*. 
We assume that we are given a nice path decomposition for G of width at most pw; we treat the 
decomposition as a sequence of bags that are ordered from left to right. To solve this problem 
we proceed as outlined above: For each partition into internal, endpoints, and unused vertices 
we take a basis for the perfect matchings on the endpoints and compute (and store) the number 
of partial solutions that are consistent with each perfect matching. We maintain and process 
this information throughout the dynamic programming; the main work is spent (unsurprisingly) 
on bags that introduce edges since this causes a rather involved recomputation of fingerprints 
(as we cannot work explicitly on separate partial solutions). 

For technical convenience our algorithm "guesses" one edge incident on a vertex of degree 
at most pw to be used in the Hamiltonian cycle. Given a nice path decomposition it can be 
easily seen that the rightmost introduce vertex bag can only introduce a vertex v of degree at 
most pw: all its neighbors must be in the current bag and no additional possible neighbors can 
be added on the right. It can be easily verified that all remaining bags, namely introduce edge 
and forget vertex bags, can be reordered freely under the constraint that no vertex is forgotten 
before all its edges were introduced. Thus, picking any edge incident on v, say {u, v}, we can 
reorder such that the last bags are: 1) introduce {u, v} (with current vertex set {u, v}, 2) forget 
vertex u, 3) forget vertex v. Our computation (from left to right) may then stop at bag 1) and 
(for some choice uj* of total weight) check the parity of the number of partial cycle covers that 
would form a Hamiltonian cycle of weight oj* when augmented with the edge {u, v}. 

Let the vertices of G = (V, E) be ordered arbitrarily, say V = {v±, ... ,v n }, and let the 
weight of any edge {vi,Vj} be given by ui(vi,Vj). We perform dynamic programming on the 
given path decomposition, proceeding from left to right (until we reach the introduce edge 
bag of the "guessed" edge). At each bag, with some vertex set B, we compute table en- 
tries t[Bo, Bi, B2, w, M] for all partitions B = Bq U B\ U B2, all integers uj £ {0, . . . , n ■ o; max }, 
and all perfect matchings M from a basis for B\ (the latter is according to Definition 3.1 with 
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a standard arbitrarily fixed ordering induced from V). Each entry contains the partity of the 
number of partial cycle covers C of the graph induced by all vertices left of and including the 
current bag and all edges introduced so far, such that 

1. C U M is a single cycle, 

2. the total weight of the edges in C is equal to u, 

3. the vertices in Bi have degree exactly i in C, 

4. and all vertices that only occur left of the current bag have degree two; we denote those 
by B t . 

We call C a (Bq, B±, B2, Bt, w)-eycle cover if it respects [2l, [3j, and |4j If it respects all four 
properties then we call it a (Bq, Bi, B2, B#, u, M)-cycle cover, i.e., if additionally the union 
with M is a single cycle. 

The main technical difficulty in the dynamic programming lies in handling the information 
stored with respect to the basis for perfect matchings of B\ , in particular when introducing a new 
edge in the path decomposition. It is crucial that we can efficiently compute a representation 
of the same information with respect to a different ordering. Intuitively, the following lemma 
allows us to change the basis of our representation. Since we apply the lemma separately for 
each partition B = Bq(jB\\JB2, set of previous vertices B%, and choice of weight uj, we state it in 
terms of a simplified table with one entry T[M] for each basis matching M . The lemma applies 
to our dynamic programming application by letting C be the set of all (Bq, B\,B^,Bi, w)-partial 
cycle covers and letting S = B%. 

Lemma 4.1. Let C denote an arbitrary set of partial cycle covers and let X denote a family 
of basis matchings on some totally ordered set S = {vq, ■ ■ ■ , wt-i} of vertices of even size. 
Furthermore, for each M £ X, let T[M] denote the number of partial cycle covers C G C 
such that M U C is a single cycle. For any X' with respect to any other ordering of S, the 
corresponding values T'[M'] for M' £ X' can be computed in time 2 t / 2 ~ l t° <yl \ 

Proof. Clearly, any permutation of the ordering of S can be achieved by at most t 2 swaps of 
two consecutive elements. Thus, it suffices to show how to move some vertex v one step to the 
"right" by swapping it with its successor in the order. We are able to show that the computation 
of any value T'[M'] requires only the contents of at most three other entries in T[-]. 

Recall the grouping of {vq, . . . ,vt-±} into vq \ v±,V2 \ ■ ■ ■ \ vts,vt-2 | «t-i. We have to 
distinguish a few cases about the position of v (odd or even) and a few of the edges in M that 
involve elements close to v in the ordering. For brevity we shorthand a little: we use a \ b,c \ 
d, e I / to denote a part of the ordering (and its groups), with a matched to b or c and / matched 
to d or e. This includes that we do not specify the positions of a and / in their groups, since 
it can be checked that they are immaterial for the discussion below (we know that one of the 
elements in the preceding and one in the subsequent group are matched like that). 

If v G {v\, V3, . . . , vt-i} (i.e., odd number and even position), then we have a | v, c | d, e | / 
(and v cannot be in the first group since that only contains vq). Consider the same ordering 
but with c and v flipped, i.e., (locally) we have a\c,v\d,e\f, and let M be a basis matching 
for that ordering. It is easy to see that M is also in the basis for the initial ordering, and 
the corresponding bitstring is obtained by inverting the bit that corresponds to the choice of 
matching a to v or c. (The bit effectively decides whether o is matched to the first or second 
element of the group \c,v\; flipping the order in the group as well as the bit cancels out.) 

The case of an odd position is, unfortunately, more involved since v moves to a different 
group. Let v be in an odd position, i.e., v S {vq, t>2, t"4, • • • , vt-2}- Thus we have a \ b,v \ d, e \ f 
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(for now let us assume that v is not the first element in the ordering). Moving v one step to 
the right results in a \ b, d \ v, e | /; let M be some basis matching for this ordering. There are 
four ways in which M can match the vertices a, b, d, v, e, / taking into account our choice of a 
and /, and they have to be treated differently (the cases are equivalent to the four matchings 
in the basis Xg for six vertices). 

i) If {a, b}, {d, v}, {e, /} E M, then M is also in the basis for the initial ordering, and with 
the same bitstring. Hence 

T'[M] = T[M]. 

ii) If {a, b}, {d, e}, {v, /} E M, then consider the following two matchings M\ and M 2 which 
are obtained by replacing the matching on d, e, v, f by either of the other two possible options 
(there are three different perfect matchings on 4 elements). 

M l := (M \ {{d, e}, {v, /}}) U {{d, v}, {e, /}, 
M 2 := (M \ {{d, e}, {v, /}}) U {{d, /}, {e, v}. 

It is easy to see that M\ and M 2 are basis matchings for the initial ordering a \ b,v \ d,e \ f, 
since all their edges connect adjacent groups. Crucially, any partial cycle cover C E C gives a 
single cycle with M if and only if this is true for exactly one of Mi and M%. To see this, it 
suffices to consider the ways in which a,b,v,d,e, f are connected (the connection through C 
can be abstracted to simply a perfect matching on these six vertices). 

T'[M] = T[M{\ + T[M 2 ] 

iii) If {a, d}, {b, v}, {e, /} E M, then we use the same edge-flipping argument, but on {a, d} 
and {b, v}. We obtain 

Mi := (M \ {{a, d}, {b, v}}) U {{a, 6}, {d, v}, 
M 2 := (M \ {{a, d}, {b, v}}) U {{a, v}, {b, d}, 

which are both in the basis with respect to the initial ordering. Thus, the desired table value 
can be computed as 

T'[M] = T[Mi} + T[M 2 ]. 

iv) Finally, we have {a, d}, {b, e}, {v, /}. In this case we use three basis matchings from the 
previous ordering, namely 

Mi := (M \ {{a, d}, {b, e},{v, /}}) U {{a, b}, {d, /}, {e, v}}, 
M 2 := (M \ {{a, d}, {b, e}, {v, /}}) U {{a, v}, {b, d}, {e, /}}, 
M 3 := (M \ {{a, d}, {b, e}, {v, /}}) U {{a, v}, {b, e}, {d, /}}. 

Again, these matchings are in the basis for the initial ordering since all edges are between 
adjacent groups. It can be checked that the desired table entry is exactly the sum modulo two 
of the table entries corresponding to these three matchings. To see this consider Table [TJ and 
identify (a, d, 6, e, v, f) with (0, 1, 2, 3, 4, 5). Then {a, d}, {b, e}, {v, /} corresponds to matching 
1, Mi to matching 6, M 2 to matchings 10 and M3 to matching 12, and the corresponding rows 
are easily seen to be linearly dependent. 

T'[M] = T[M{\ + T[M 2 ] + T[M 3 ] 

Since there are 2 t//2_1 perfect matchings in the basis for any ordering of S and each entry 
takes only 0(1) operations the claimed total time follows by computing all intermediate tables 
(only the most recent one needs to be stored). □ 
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Now we can return to the description of our algorithm. The dynamic programming on 
the path decomposition proceeds from "left" to "right", using the table of the previous bag, 
denoted t'[], to compute the table of the current bag, denoted t[]. We use t' uv and t uv for the 
corresponding tables that are obtained by changing ordering and basis in such a way that u 
and v are the last two vertices in the total order. For introduce and forget vertex bags there 
is not much work required since by themselves vertices do not affect our partial solutions. The 
main work lies in the computations required for the introduce edge bags: Partial solutions that 
use the new edge have a different set B\ of degree- 1 vertices, which comes with a different basis. 

First bag: For the first bag the vertex set B is empty and we only get a single trivial table 
entry 

£[0,0,0,0,0] = 0. 



Introduce vertex bag: We have a current bag with vertex set B that introduces a vertex v. 
Accordingly the previous bag has vertex set B \ {v}. Furthermore, the set of vertices that occur 
only on the left is the same for both bags. Clearly, in the subgraph given by all vertices of the 
current and preceding bags plus all edges introduced so far no partial cycle cover can include 
edges incident on v. Thus 

t[Bo, Bi, B 2 , •, •] = 0, 

for all partitions B = Bq U B\ U B 2 with v G B\ U B 2 . For partitions with v G Bq and any choice 
of weight to and i?i-basis matching M we already have the correct number in the table for the 
previous bag, namely 

t[B , B!,B 2 ,lo, M] = t'[B \ {v}, Bx,B 2 ,u, M]. 



Forget vertex bag: We have a current bag with vertex set B that "forgets" vertex v ; the 
previous bag has vertex set B U {v}. Let Bi be the vertices that occur only left of the current 
bag and note that Be \ {v} occur only left of the previous bag. Fix any partition B = Bq U 
B\U B 2 , a matching M from the .Bi-basis (with standard ordering), and a weight to. If C is 
a (So, Bi,B 2 , Bi, to, M)-cycle cover, then v must have degree two in C by definition. Conversely, 
any (Bq, Bi, B 2 U {v}, B^ \ {v}, to, M)-cycle cover C (whose number modulo two is stored in the 
previous bag) is also a (Bq, B\, B 2 , Bi, to, M)-cycle cover. Hence 

t[B , B u B 2 ,u, M] = t'[B , B U B 2 U {v},u, M}. 



Introduce edge bag: We have a current bag with vertex set B that introduces an edge {u, v} 
with u, v G B. The previous bag has the same vertex set B (but its partial solutions cannot 
make use of {u, v}), and both bags have the same set Bi of vertices that occur only left of them. 
First, we compute from the table of the previous bag the table t' uv []; i.e., the information 
is represented with respect to a basis that has u and v as the last two vertices in the ordering 
using Lemma 4.1 We will proceed by computing t uv \\ which will subsequently be transformed 
to t[}; this completes the procedure for an introduce edge bag. 

We explain the computation of t uv [BQ, B\, B 2 ,to, M] from t' uv \} for an arbitrary partition B = 
BqL) B\L) B 2 , integer weight to, and perfect matching M; the latter is from the basis for B\ with 
modified ordering of V (u and v are last two elements). We need to consider different cases 
depending on which sets Bi contains u and v. 
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i) u G Bq or v G Bq: In this simple case no (Bq, B\, B 2 , w)-cycle cover may contain the 
edge {u, v} and the same is true for all perfect matchings M on B\ (in particular for those from 
a basis). Thus, directly from the definition of our tables we get that 

t uv [B , B 1} B2,u,M]= t' uv [B , Bi,B 2 ,uj,M], 

for all basis matchings M. 

ii) u, v G B\: Since M is a perfect matching from a basis on B\ it may contain the 
edge {u, v}. If {u,v} G M then no (Bq, B\, B2, Bg, ui, M)-cycle cover C can contain {u, v} 
(since the two copies of {u, v} would already give a separate cycle). (In the final introduce edge 
bag we do not need to perform this computation since we only check for a solution.) We already 
know the parity of the number of (Bq, B±, B2, Bg, u, M)-cycle covers avoiding the edge {u, v} 
from t' uv [] and get 

t uv [B , B ± ,B 2 ,cu, M] = t' uv [B , Bi, B 2 ,oj, M]. 

Otherwise we have {u,v} £ M which implies that {p, u},{q,v} G M for some p,q G B\\ 
{u, v} (with p / q) since M is a perfect matching of B±. As before we already know the parity 
of the number of (Bq, B\, B 2 , Bg, 10, M)-cycle covers that do not contain {u, v} from t' uv []. 

Consider a (Bq, B\,B 2 , Bg, cj)-cycle cover C that does contain {u, v}. Thus CUM contains a 
cycle S = (xi, . . . , x r ,p, u, v, q, x\). Intuitively, we want to compare this to partial cycle covers 
without {u,v}; we effectively move the edge into the matching: Let C := C \ {{u, v}} and 
let M' := (M \ {{p, u}, {q, v}}) U {{p, q}} (the additional edge in M is modeled by contract- 
ing p, u, v, q into just p, q). Observe that C is a (BqU{u, v}, B±\{u, v }, B 2 , Bg,uj — u(u, v))-cycle 
cover. Clearly, C U M' contains a cycle S' = (x\, . . . ,x r ,p,q,x±), and all further cycles (if they 
exist) are the same as in C U M. Thus C is a (Bq, B±, B 2 , Bg, u, M)-cycle cover if and only 
if C is a (Bq U {u, v},B\\ {u, v},B 2 , Bg, u — oj(u, v), M')-cycle cover. This fact is useful for the 
computation, provided that M' is in the basis for B\ with modified ordering, which will check 
next. 

Let us see that M' is indeed in the basis for B\ \ First of all, since u and v are at the end of 
the ordering of B\ , when we remove these two vertices the remaining vertices keep their current 
ordering. Second, due to the structure of bases for B±, we see that the edges {p, u}, {q, v } G M 
imply that the ordering of B\ ends either with . . . ,p,q,u,v or with . . . ,p,r,q,u,v (to see this, 
recall the grouping of the elements and the fact that the basis matchings are exactly all ways 
of pairing up elements from adjacent groups). Accordingly the ordering for B\ \ {u,v} ends 
in . . . ,p, q or . . . ,p, r, q. If we have . . . ,p,q,u,v then M = X(\Bi\, alO) and M' = X(\Bi\ — 
2,al), for some bitstring a. If we have . . . ,p,r,q,u,v then M = X(\B\\,a§Q) and M' = 
X(\B\ \ — 2, aO). In both cases, M' is part of the basis for B\ \ {u, v} (in fact, its corresponding 
bitstring is the same as the one for M minus the last position). Thus, taking into account the 
contribution for cycle covers without {u, v}, we can compute t uv [BQ, B\, B 2 , oj, M] as follows 

t uv [Bq , B l , B 2 , uj , M] = t' uv [Bq,B u B 2 ,u, M] 

+ t' uv [B U {u, v}, B x \ {u, v}, B 2 ,uj - u(u, v), M'\. 

iii) u G B\ and v G B 2 : We get that {u,v} ^ M, but there must be some p G 
B\ \ {u} with {p, u} G M since M is a perfect matching of B\. Again we get the contri- 
bution of t' uv [BQ, B\,B 2 ,oj, M] for all (Bq, B\,B 2 ,Bg, u>, M)-cycle covers that do not contain the 
edge {u, v}. 

If a (Bq, Bi, B 2 , Bg, w)-cycle cover C contains the edge {u, v}, then CUM contains a cy- 
cle (x\, . . . , x r ,v, u,p, xi). Again, we effectively move the edge {u, v} "into" M: Consider C' := 
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C \ {{u, v}} and M' := (M \ {{p, u}}) U {{p, v}}. Observe that C is a (B U {u}, (B x \ {u}) U 
{v}, B2\{v}, Bi, uj—oj(u, -u))-cycle cover. Furthermore, C'UM' contains a cycle (xi, . . . ,x r ,v,p,xi), 
and all further cycles (if there are any) are the same as in CUM, Thus M U C is a single 
cycle if and only if M 1 U C is a single cycle. Again we need to check that M 1 is in the basis 
for (B\ \ {u}) U {v}. In the present case this is straightforward: Since u and v are last in 
the modified ordering of V, they both occupy the last position on B\ and (B x \ {u}) U {v} 
respectively. Since M and M' differ exactly by replacing u with v, we get that M 1 is in the 
basis for (B\ \ {u}) U {v} and, in fact, corresponds to the same bitstring as M for B\. We get 

t uv [B , B u B 2 , u, M] = t' uv [B , B h B 2 ,u, M] 

+ 4, [Bo U {u}, {B l \ {u}) U {v}, B 2 \ {v},oj- u(u, v),M'}. 

iv) u G B 2 and v € B\: This case is symmetric to the previous one. Despite u being second 
to last and v being last in the modified ordering, this requires no change of argumentation since 
we always have only one of them in the set of matched vertices. 



v) u,v G B 2 : Clearly no perfect matching M on B\ can contain {u,v} or any other edge 
incident with u or v. However, (Bq, B\, B 2 , Bi, w)-cycle covers C may contain the edge {u,v}. 
The parity of the number of such partial cycle covers that do not use {u, v} and that are 
consistent with M is already stored in t' uv [Bo, Bi, B 2 , u, M]. 

Let us consider a (Bq, B\, B 2 , Bi, w)-cycle cover C that does contain {u,v}. Let C := C \ 
{{u,v}} andM' = MU{{u,v}}. Observe that C is a (Bq, BiU{u, v}, B 2 \{u, v}, B e , uj-uj(u, v))- 
cycle cover. It is easy to see that C U M is a single cycle if and only if C U M' is a single cycle, 
since CUM = C'UM'. 

Let us check that M 1 is in the basis for B\ U {u, v} with modified ordering. Let the end 
of the ordering of B\ be ...,p,q,r. If M = X(\Bi\,aO) then {p, r} G M, and we get M' = 
X(|Bi| +2,o01). Else, if M = X(\B x \,al) then {q,r} G M and we get M' = X(\B X \ +2, all). 
(Note that we basically only need to append a 1 to the bitstring for M to enforce that p or q 
is matched to r and not to u; the edge {u, v} G M' is then implied.) Thus M' is indeed in the 
basis for B\ U {u, v} with modified ordering. We get 

t uv [B , Bi,B 2 ,u, M] = t' uv [B ,B 1 ,B 2 ,co, M] 

+ t' uv [B , B x U {u, v}, B 2 \ {u, v},u - oj(u, v),M U {{u, v}}}. 

Now we have handled all cases for computing t uv [Bo, B\, B 2 ,u, M]. It remains to apply 



Lemma 4.1 again to compute the table t\\ that represents the data with respect to the standard 



ordering. This completes the necessary work for an introduce edge bag. 



Runtime. Let us now analyze the time required to compute all table entries for one bag. 
We already know how to compute t[Bo, B x , B 2 ,uj, ■] in time 2l Bl l //2_1 |f?i| c '( 1 ), spending effec- 
tively l-Bil^ 1 ) per table entry (each entry corresponds to one basis matching on B\). For every 
weight 1 < u) < nw max the number of tuples (Bq, B x , B 2 ,M) where Bq,B\, B 2 partitions of the 
bag B and M is one of the 2^ Bl ^ 2 ~ 1 basis matchings equals (2 + \/2)l B l by the multinomial 
theorem, and hence we have proved the following: 

Lemma 4.2. There is an algorithm that given a graph G along with a path decomposition of 
width pw computes the table entries t\\ corresponding to the rightmost introduce edge bag of the 
path decomposition in (2 + V2) w w majX (n ■ pw)^ 1 ) time. 

Now, we can wrap up by completing our algorithm. 
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Theorem 4.3. There exists a Monte-Carlo algorithm that given a graph G along with a path 
decomposition of width pw solves the Hamiltonian cycle problem in (2 + \^2) p ^(n ■ pw)^ 1 ) time. 
The algorithm cannot give false positives and may give false negatives with probability at most 
1/2. 

Proof. The algorithm is the following: First guess an edge {u, v} in the solution that is inci- 
dent on the vertex v introduced in the rightmost introduce vertex bag, and reorder the path 
decomposition as described above. (The guessing is of course implemented by trying all of 
these at most pw many edges.) Assign for each edge e G E a weight uj{u) G {1, 2, . . . , 2\E\} 



uniformly and independently at random. Then use Lemma 4.2 to compute the table entry 
t[0, {u, v}, 0, {{u, v}}, w] corresponding to the rightmost introduce edge bag for every weight 
w < n • u> max < n2\E\ in 2 t / 2 (nt) oi ^] Return YES if one of these values is 1 and NO otherwise. 

This procedure clearly runs in the claimed running time. If the algorithm return YES, 
there exists a Hamiltonian cycle including the edge {u, v} since there exists a Hamiltonian cycle 
by the definition of t. Otherwise, suppose that there exists a Hamiltonian cycle. Then with 



probability at least 1 — |i£|/2|.E7| = 1/2 we have by Lemma 3.9 that u isolates the family of all 
Hamiltonian cycles of G. Hence with probability at least 1/2 there exists a unique minimum 
weight Hamiltonian cycle H minimizing uj(H). Let {u, v} be the edge incident to v contained 
in H. Then we have that t[^,{u,v},$,{(u,v)},u(H) — u({u,v})] = 1 since H is the only 
Hamiltonian cycle of weight to(H). □ 

Further results 



Combining the ideas from the proof of Theorem 4.3 with Lemma 4.4 and a trick from |15j . 
we will now give an application to fast exponential time algorithms for the Hamiltonian cycle 
problems on graphs of degree at most 3. To obtain this we require the following result: 

Theorem 4.4 ([IH]). For any e > there exists an integer n t such that for any graph G with 
n > n t vertices, 

rr,\ ^ 1 1 13 
pw(G) < -n 3 + -n 4 + —n 5 + n> 6 + en, 

b 3 60 

where n{ is the number of vertices of degree i in G for any i E {3, . . . , 5} and n>6 is the number 
of vertices of degree at least 6. 

This theorem is constructive, and the corresponding path decomposition (and, consequently, 
tree decomposition) can be found in polynomial time. The observation from |15j is that on cubic 
graphs, we know that the degree of a vertex in a partial cycle cover of the type used in the 
our table entries cannot be 2 if at most 1 incident edges is introduced, and it cannot be if at 
least 2 edges are introduced. Hence these states can be safely ignored since they will be empty. 
Summarizing, we obtain the following: 

Corollary 4.5. There exists a Monte-Carlo algorithm that given a graph G of maximum degree 
three solves the Hamiltonian cycle problem in (1 + \^2) n ^n c '^ < 1.1583 n n°( 1 ' time. The 
algorithm cannot give false positives and may give false negatives with probability at most 1/2. 



Again the proof idea similar to proof of Theorem 4.3 we can also rather easily obtained the 
following result: 

Theorem 4.6. There exist algorithms that given a graph G and integer k finds a path of length 
k in time n(2 + \^2) pv (k + pw)^ 1 -* time if a path decomposition of width pw of G is given. 
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The idea is to slightly modify the table entries defined in this section: In the definition 
we considered partial cycle covers visiting all already forgotten vertices. We allow to not visit 
some vertices and keep track of the number of visited vertices with an additional counter, and 
hence the formula in the forget bag has to be altered to t[Bo, B\, B2, u, M, £] = t'[Bo, B\, B2 U 
{v},oj,M, £] + t'[Bo U {v}, Bx, B%, to, M, £] where £ is the additional counter. The additional 
counted can be bounded by the observation that a non-zero entry with £ > pw • k directly 
guarantees a solution. 



5 Lower bound for Hamiltonian cycle 

In this section we prove that the runtime of our algorithm from Section [4] is essentially tight. 
We show that no algorithm can achieve a significantly better dependence on the pathwidth of 
the input graph (even at the cost of a larger polynomial factor in the input size) , without giving 
also a breakthrough result for solving CNF-Sat; this is expressed by the main theorem of this 
section. 

Theorem 5.1. If any algorithm solves the Hamiltonian cycle problem in (2 + y/2 — e)P w | V"^ 1 ) 
time, then there exists an algorithm that solves the CNF-Sat problem in (2 — e) n mP^' time. 

Recall the definition of the Matching Connectivity matrix Tit- All rows and columns are 
indexed by perfect matchings M\ , M2 of the complete graph on t vertices and Tit [Mi , M2] = 
[M1UM2 is a Hamiltonian cycle]. Our proof, namely a reduction from CNF-Sat, uses the lower 
bound on the rank of Tit, in particular, the fact that the family of basis matchings Xj from 
Section [3] induce a permutation submatrix in Tit (see Proposition 3.3). 



5.1 Gadgets 

In this section we will introduce three gadgets used for the final construction. Since our gadgets 
accept some parameters we will call a graph that contains some gadgets an annotated graph. 



Vertices with labeled incident edges. The following two gadgets allow to label incident 
edges of a vertex v with a labeling function X v , ensuring that every Hamiltonian cycle enters 
and leaves a vertex with edges of the same label. 

Definition 5.2. A label gadget is a pair (y, X v ) where A^ is a labeling of the edges incident to 
v. A Hamiltonian cycle C is consistent with a label gadget (v,X v ) if X v (e) = X v (e') where e,e' 
are the two edges of C incident to v. 

The first gadget is for a vertex whose incident edges have only two distinct labels, denoted 
by dotted and normal lines in figures. It is shown at the left-hand size in Figure [2} We have 
to ensure that any Hamiltonian cycle contains exactly two edges of the set of edges leaving 
the gadget and that they are of the same label, and this can be seen to hold by a simple case 
analysis: If the cycle enters v\, it must continue with V2,v%. Then it cannot leave the gadget, 
because then it is impossible to visit all six remaining vertices. Hence it must continue with vq 
and then v&, V4, vg, vs, v-j is forced. The cases where it enters at a different vertex are symmetric. 

Let us continue with the multilabel gadget, which we will first formally describe. Given a 
vertex v in a graph G with incident edges X and a labeling A^ : X — >■ [k], we obtain G' by 
replacing v as follows: 

• Create a cycle of 3/c vertices consecutively denoted by vf, v\, vf, v%, v\, v%, ■ ■ ■ , 
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Figure 2: Gadget for vertices with labeled incident edges. On the left-hand side the 2-label 
gadget is shown. The two label classes are denoted by normal and dotted lines. On the right- 
hand side the multi-labeled gadget for six labels is shown. The black vertices represent 2-label 
gadgets and the labels of their incident edges are denoted by normal and dotted lines. A number 
i next to a leaving edge indicates the edge represents an edge in the original graph with label i. 

• For every 1 < i < k, and every vertex w / v incident to an edge e G X such that A„(e) = i, 
add the edges vfw and v\w. 

• For every 1 < i < k add a 2-label gadget consisting of vertex gi, edges givf, giv\ labeled 
by 1 and the edges g%vf, l5 giV b i+l labeled by 2. The gi are called guard vertices. 

An instructive example with k = 6 is shown on the right-hand side of Figure [2} 

Lemma 5.3. There is a Hamiltonian cycle of G consistent with (v, X v ) if and only if G' has a 
Hamiltonian cycle. 

Proof. Note that in a Hamiltonian cycle a guard vertex gi is either adjacent to both vf and v\ 
or to both vf +l and v!f +1 , and all the pairs can be adjacent to at most one guard vertex since 
otherwise there is a cycle on four vertices. Moreover, since each vertex vf is of degree two, any 
Hamiltonian cycle of G' contains the edges v\vf for 1 < i < k, as well as the edges vfvf +1 for 
1 < i < k and finally the edge v?vf. Therefore, since all guard vertices must be visited, any 
Hamiltonian cycle of G' contains only two edges leaving the gadget and they must be incident 
to vf and v\ for some i. 

For the forward direction, let there be a Hamiltonian cycle of G that enters and leaves at 
edges incident of v with label i. Then we have to show there is a Hamiltonian path P in the 
gadget from vf to v\, ignoring edges leaving it. Notice that such a P exist in which (i) gj is 
adjacent to and v b for j < i, and (ii) gj is adjacent to and Vj +l for j > i. □ 

Induced subgraph gadget. We will now discuss the induced subgraph gadget, which is 
shown in Figure [3j The function of the gadget is described as follows: 
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Figure 3: Example of the induced subgraph gadget, with X = {t>i, . . . , uio} and F = 
{{viv 3 , v 2 vq, viVq, v 4 v 5 }, {v 3 vq, v 5 v 7 , v 6 v 8 , v 9 v w } , {vgv 10 , v 8 v 9 , v 7 v w }} . The double triangles at 
the Ijr indicate that in fact there are edges to every vertex in I jr. Label 1 edges are drawn 
dotted while label 2 edges are drawn normally. 
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Definition 5.4. An induced subgraph gadget in a graph G = (V, E) is described by a tuple 
{X, a, b, F) where X CV, a,b eV\X, F C 2 £ ( x > x ) and F. A Hamiltonian cycle G of G 
is consistent with (X, a, b, F) if ab £ C and C n E(X, X) £ J 7 . 

The gadget is implemented with the following construction, obtaining a modified graph G' 
from G: 

• remove all edges E(X, X) from the graph G, 

• add a set /iw, i.e. an independent set of jJ 7 ) vertices to G, and make all its vertices 
adjacent to both a and b, 

• Let J 7 = {Fi, . . . , i^}, and for % = 1, . . . , £ do the following: 

- Let Fi = {ei, . . . , e z }. 

- Add a path of 2-label gadgets Pi = {pj, . . . ,p\ }, with all edges having label 1. 

— For every vertex v of the independent set i^, add an edge with label 1 to p\ and 

p\ ■ 

— For j = 1, . . . , add edges p^x and p^y with label 2, where ej = xy. 

We will now prove the function of the gadget. The vertices a and b are not useful for using 
the gadget but are required to be able to implement it. 

Lemma 5.5. There exists a Hamiltonian cycle C in G that is consistent with (X, a, b, F) if and 
only if there exists a Hamiltonian cycle in G' . 

The intuition of this is that in a Hamiltonian cycle of the graph the part of the cycle that 
goes from a to b visits all vertices from the independent set and all paths except the one 
corresponding to corresponding to the element of F that is chosen. 

Proof. For the forward direction, let Fi 6 F be the intersection of G with E(X,X). Then the 
Hamiltonian cycle in G can be extended to a Hamiltonian cycle in G' by replacing every edge 
ej £ Fi with two edges up? and p\v where ej = uv and by replacing the edge ab by the path 
a, 9i, Pi, 92, P2, ■ ■ .,Pi-i,9i,Pi+i,gi + i, ■ ■ .,P\j:\g\j:\,b, where g u . . .,g\jr\ is an arbitrary ordering 
of 

For the reverse direction, let C' be a Hamiltonian cycle of G' . Observe that since is 
independent, the Hamiltonian cycle C' contains exactly 2\F\ edges between and Ng'(I\f\) = 

{pj,Pi ■ 1 < i < U {a, b}. By the 2-label gadgets, we know that if for some 1 < i < \F\ 
the Hamiltonian cycle C' contains an edge between 1^ and {p},p^}, then C' contains exactly 
one edge between Ljr\ and p\ as well as exactly one edge between Iyp\ and p\ Fl ^- However C' 

cannot contain 2\F\ edges between Ljr\ and {p\,p\ : 1 < i < \F\}, because then it would 
not visit any other vertices. Furthermore C' cannot contain two edges between a (or b) and 
-Tiw, since then it would have to contain either zero or two edges between b (or a) and Iiji, 
making it impossible to visit the original vertices of G. Consequently G' contains exactly one 
edge between a and exactly one edge between b and /iw, and exactly two edges between 
Pi and for \F\ — 1 indices i G [\F\]. So there is exactly one index 1 < io < \F\ such that 

p\ and are not connected with an edge of Since C' visits all the vertices of Pi by 
label 2 edges only, we can obtain a Hamiltonian cycle G in G by removing all the edges of the 
gadget and adding the edge ab together with edges of Fi . Note that the Hamiltonian cycle G 
is consistent with (X, a, b, F) as it does not contain any edge of E(X,X) \ Fi since the graph 
G' has the edges E(X, X) removed. □ 
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Figure 4: The total construction for the Hamiltonian path lower bound. The Kg denote cliques 



of f3 vertices. The black circles Cj are multilabel gadgets (see Figure 
are induced subgraph gadgets, used as further illustrated in Figure [5 



and the rectangles Gi 



5.2 Construction 

We will now describe the construction of the reduction, see also Figure |4j After the formal 
definition we will provide some intuition. Our construction is parameterized by two integer 
constants j3 and 7, which will be chosen later, such that 



V2 1 > 2 7 . 



(1) 



Assume we are given a CNF-formula = C\ A . . . A C m on variables x±, . . . ,x n with n being 
a multiple of 7; let q = ^. Partition the set {xi, . . . , x n } into n/7 blocks of size 7, denoted 
Xi, . . . , X n /^. Also, denote Xi = {xn, . . . , Xj 7 }. Intuitively, we will represent the 2 7 assignments 
of a block of variables by the states of group of vertices in a bag of the to be constructed path 



decomposition. In Subsection 5.5, we will discuss how to fix f3 and 7 such that the number of 



these states, represented on the left hand side of Equation ([!]) is at least 2 7 . 

In the following we will use the family of matchings defined in Definition 3.1 For ease 
of notation we denote X(S') for the family obtained from X|5| by replacing the elements of Ut 
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Figure 5: An edgeset from T in an induced subgraph gadget G\\ is shown here. For convenience, 
Vi^o and ^2,10 are denoted by v\q and V2o- The edgeset shown is constructed from the coloring 
s with s(vii2) = s(vm) = 0, 5(^114) = s(viis) = s(vhq) = s(vns) = 2 and the matching 
(^111,^113), (^119,^11,10)- 

one to one with elements of S (using an arbitrary but fixed ordering), for any set S of even 
cardinality. 

1. For every i = 1, . . . , q, j = 1, . . . , m + 1 and k = 1, . . . , (3 add vertices for i = 1, . . . , q 
and j = 1, . . . , m + 1 denote Vij = {%i, . . . , %^}. 

2. For every i = 1, . . . , q add a clique on /3 vertices K%, and add an edge between all vertices 
of the clique Kp and all vertices from V% t i. 

3. For every i = 1, . . . , q add a clique on f3 vertices K%, and add an edge between all vertices 
of the clique Kl and all vertices from Vi )m +i. 

4. For every i = 1, . . . , q — 1 add a vertex and make it adjacent to all vertices of and 
K p + . Furthermore, a vertex y' q is added and made adjacent all vertices of K q ^ and Kp. 

5. For every i = 1, . . . , q and j = 1, . . . , m we use an induced subgraph gadget Gij (the special 
vertices a*- and 6*- from Definition 



5.4 



are included in Gij) as follows (see also Figure pi): 
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(a) Add four vertices a*-, c*-, dj and the edge (a*-, &*•), where the edge (a*-, 6*) is guaran- 
teed to be in any Hamiltonian cycle (by for example subdividing it) as required by 
the gadget. 

(b) The set X consists of the vertices Vij U and the vertices c*-,d*. 

(c) Given a vector s G {0,1,2}" and a matching M G X(s _1 (l)) we denote suj\ as 
the coloring of Vij and M^j^ as the matching of Vij obtained by identifying the 
sets {1, . . . ,/3} and Vij. (Equivalently, we may take suj\ G {0, 1,2}^ and Muj\ G 

x(^ } (i)).) 

(d) For a every vector s G {0,1,2}^ and a matching M G X(s _1 (l)), add an edgeset 
r]ij(s, M) constructed as follows to the family T used for the induced subgraph gadget: 

i. add the edges of the path ai,...,ai where {ot\, . . . ,ae} are the elements of 
s7\(2) sorted ascendantly (the ordering is in fact immaterial so long as we have 
a simple path through these vertices); 

ii. add the edges of the path /?i , . . . , where . . . , fy/} are the elements of 
s7^- +1 \(0) sorted ascendantly; 

hi. let e = (yi, 2/2) be the edge of incident with the smallest element of s~zK (1); 

iv. add all edges of Muj\ not equal to e; 

v. let M 1 G X(s _1 (l)) be the unique matching such that M U M 1 is a Hamiltonian 
cycle (this exists and is unique due to Proposition 3.3); 

vi. let e' = ( zi, z 2 ) be the edge of incident with the smallest element of 
s7.\(l)] 

vii. add all edges of . +1 v not equal to e'; 

viii. add the edges (y 1 ,a t j ),(y 2 ,c l j ),(c l j ,a 1 ); (z x , 6J-), (z 2 , d l j), (d},/?i); (ai,Pe). 
6. Fix arbitrarily an injective mapping 

Tp : {0, l} 7 -»• {(s, M) s G {0, 1, 2} 13 A |s _1 (l)| is even A M G X(s" 1 (l))} . 



7. For j = 1, . . . , m: 

(a) Add a multilabel gadget Cj. 

(b) For every i = 1, ...,q and every partial assignment x G {0,1}^* of the variables 
Xi that satisfies Cj, consider the first vertex p of the path in the induced subgraph 
gadget Gij corresponding with the element T)ij(ip(x)) of J 7 . 

(c) Let (p, a) be an edge with label 1 in the construction of the induced subgraph gadget. 
Add an edge (p,Cj) and (cj,a) such that with respect to p and a the edges again 
have label 1 and with respect to the multilabel gadget Cj the edges have a label that 
has not been used yet in the interval [j ■ T 1 + 1, (j + 1) • 2 7 ], which must exist since 
there are at most this many assignments of Xi. 

The intuition behind the construction is as follows: Vij are blocks of vertices whose joint 
states encode a joint assignment of X;. In Step 2 — 4 vertices are added to ensure that a set 
of disjoint paths visiting all vertices in the remainder of the graph can be completed into a 
Hamiltonian cycle. In Step 5 we create graphs Gij on Vij and Vij+i such that any for any set of 
these disjoint paths, the intersection of Gij and Gij+\ is similar. To do this we use the induced 
subgraph gadget and for every state (s, M) we allow exactly one edgeset in T that induces state 
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(s, M). In Step 6 we fix an encoding tp of partial assignments in states of the group of vertices 
Vij. As mentioned before, this encoding exists due to Q. Finally, we check in Step 7 whether 
clause j is satisfied by adding a multilabel gadget Cj that can be visited if and only if in some 
graph dj the intersection X from the family T induces a state (c, M) = rj^ (X) such that the 
partial assignment ^^(s^M) satisfies clause Cj. The labels to the clause gadget are chosen 
consecutively in 7c to ensure adding the labels only increases the pathwidth by a constant. 

5.3 From a satisfying assignment to a Hamiltonian cycle 

Suppose that x £ {0, l} n satisfies the formula (p. Then we first claim that for every i = 1, . . . , q 
the set 

m 

S=\Jr Hj (i>(X i )) (2) 

is a set of disjoint paths in the constructed graph G, in which all vertices of UjLiGij are visited 
and all endpoints are in Vii U Vi jm +i. To this end we first need to zoom in at the application 
of the induced subgraph gadget in Step 6. of the construction (see also Figure [5]). Note that 
effectively, in a set T]ij(s,M), the first vertex y\ of s/^-\(— 1) is connected to the first vertex z\ 
of S(jj +1 }(— 1) through Recall that y<i is the second vertex of the lexicographically first 

edge of M^jy Suppose that t/2 = vijk- 

Then, to show the claim on £, it is sufficient to show that the second vertex yi is connected 
to ft+ij by a path with internal vertices s^(2) U s^*^ ({0, 1}) U c*,(i*-. This follows by 
construction: in Vi+i we use two matchings on the vertices with label one that are guaranteed 
to form a Hamiltonian cycle so since we left out the edges (z\, Z2) in one matching and the edge 
(^i,«ij+i,fe) they must now form a path between the z^ and Vij+i^- It is easily seen that the 
further extension of the path meets the claim. 

Now we alter the set £ to a set £' to visit all clause multilabel gadgets as follows: since 
x satisfies the formula, we know by construction that for every j = 1, . . . ,m there exists i = 
1, . . . , q such a variable of Xi is satisfying Cj. We take one such i and we know that in the path 
created in the induced subgraph gadget corresponding to the element r](ip(Xij)) of T we added 
in Step 7.C the two edges to the clauses gadget with the same label between a vertex p of X 
used and the first vertex a of the created path. Since we know (p, a) £ £ we can safely replace 
0,<t) with (p,Cj) (cj,a). 

Hence now £' is a set of disjoint paths, with endpoints in U^V^i U Vi jm+ i, all vertices 
ci, . . . , Cm and vertices from UijGij are visited, and by construction it is consistent with all the 
multilabel and induced subgraph gadgets. After contracting the internal points of these paths, 
it is easy to see we are left with a graph that has a Hamiltonian cycle: for every i = 1, . . . , q 
we can use the cliques to extend the paths from £' to a path between two vertices of K% that 

contains all paths from £' with endpoints in uf^V^i U Vi t m+i that visits U Kl and does 
not visit y[_i and y[ such that the vertices y[ then can be used to connect these paths into a 
complete Hamiltonian cycle. By construction it is easily seen to be consistent with all gadgets. 

5.4 From a Hamiltonian cycle to a satisfying assignment 

Suppose we are given a Hamiltonian cycle C of G, consistent with all the induced subgraph and 
multilabel gadgets. 

Contract edges to clause gadgets obtaining C' . By the multilabel gadget and the labeling, 
the edges resulting from the contraction will be contained in an induced subgraph gadget. Let 
us denote Xij for the vertex set of Gij minus etjj , bij . 
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Lemma 5.6. For every i and j it holds that 



rj£{Cf n {X t3 x X l3 )) = r£ +1 (C n (X w x X iJ+1 )). 

Proof. Let (s, M) = ^(C n (X tj x Xy)) and (s',M') = 77^ (C" n (%i x Xy+i)). First, 
we know that for every Sj + = 2 since otherwise the vertices of V{ will not have degree 2 in 



C". Hence, we know that = (s') _1 (l). If M 7^ M' then we know by Proposition 3.3 that 

M U M' is not a Hamiltonian cycle since the matchings are by construction from X(s _1 (l)). 
Hence we have that M U M' contains two subcycles and the one avoiding the element 1 will 
give a subcycle of C", contradicting that C is a Hamiltonian cycle. □ 

Hence, the cycle C must represent an assignment of the variables x±, . . . ,x n . Furthermore, 
since C visits all multilabel gadgets, it must represent a satisfying assignment by construction. 

5.5 Pathwidth and efficiency bound 

We will first argue that the pathwidth of the constructed graph (after expanding all gadgets) 
G is " (3 + /(A 7) for some function / and then set the parameters f3 and 7 in order to prove 
the theorem. 

Let k = l,...,q and j = 1, . . . , m. In the multilabel gadget Cj, recall there is a vertex vf, 
adjacent to some vertex in the induced subgraph gadget G\j. Similarly, let t be the maximum 
integer such that in this multilabel gadget, is adjacent to a vertex in one of the induced 
subgraph gadgets G\j , . . . , G K j . Then note that 

k q 
S.,j = ((jV l3 )U( U V itj+1 )U{vi,g t ,vt}, 

8=1 n=i+l 

is a separator since the labels in the multilabel gadget are ordered such that labels coming from 
the induced subgraph gadget are adjacent. Then the claimed upper bound on the pathwidth 
follows easily by using the above separators as bags consecutively for every j = 1, . . . , m and 
K = 1, . . . , q since the size of every union of two consecutive separators is ™j f3 + f(f3, 7) and 
the components to the left of £1,1 and to the right of S q , m are of constant pathwidth. 

Now we will discuss how to fix the constants (3 and 7. Note that by splitting summations 
we obtain 



E (• ? • )v^ x = E ^2 

«0+H+»2=P 0<ii</3 
ii is even iieven 

and since for every odd i±, we have y / 2* 1 2^ - * 1 < \/2 ll 2^~ n ~ 1 we know by the multinomial 
theorem that 

il is even 

Hence, we satisfy Equation jlj by setting /3 = ^2+^) ^ ^' ^ unnm § the assumed algorithm that 
solves the Hamiltonian cycle problem in (2 + \[2 — e) pw | Vp^ 1 ) algorithm then takes time 



(2 + V2- e )M( 1 ^ +1 ) <2 lg(2+ ^- e) ^( 1 ^ +1 ) =2 (n+1) ( 



lg(2+y^2- £ ) lg(2+V^-e) ' 

ig(2+\/2) -y 



So choosing 7 large enough such that ( ^7^75^ + lg ^ 2+v ^ e ^ ) < 1 is sufficient and gives a 



lg(2+%/2) ^ 7 



(2 — ^"m ' 1 ' algorithm for CNF-Sat, concluding the proof of Theorem 5.1 
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6 Conclusions 



We have presented a set of matchings X n , which forms a basis of the Matching Connectivity 



matrix H n . In particular we have obtained a factorization theorem (Theorem 3.4) which shows 
an explicit way of expressing any perfect matching as a linear combination of matchings of X n . 
As a consequence we obtained deterministic algorithms for computing the parity of the number 
of Hamiltonian cycles in undirected graphs and directed bipartite graphs in 0(1.888"), which 
together with the Isolation Lemma lead to Monte Carlo algorithms solving the decision versions 
of Hamiltonicity within the same running time. Moreover, using the basis X n , we presented 
an algorithm which given an undirected graph on n vertices along with a path decomposition 
of width at most pw, decides Hamiltonicity in (2 + v / 2) pw ra c '( 1 ) time. Somewhat surprisingly we 
use the same tool, i.e. the basis X n , to show by an involved reduction from CNF-Sat that our 
bounded pathwidth algorithm is optimal under the Strong Exponential Time Hypothesis. 

Our results lead to several natural open problems. Can the basis X n be used to obtain a 
deterministic 0((2 — e) n ) time algorithm for Hamiltonicity? Can we handle directed graphs 
without the bipartiteness assumption? Can we extend our bounded pathwidth algorithm to a 
bounded treewidth algorithm with the same complexity? 

Finally we would like to note that the row space of the Matching Connectivity matrix T~L n 
clearly has several different bases. We have investigated a particular one, which proved to have 
several interesting properties and applications, however there might be different ones which are 
also worth exploring. 
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A Omitted proofs of Section [3] 
A.l Proof of Theorem 13.41 



This section provides a proof for Theorem 3.4, which itself implies that each family Xj gives a 
basis for the Matching Connectivity matrix T-Lf The proof will be by induction. As a first tool we 
define a projection that maps any perfect matching of Ut to a certain perfect matching of Ut- 2 ; 
this is controlled by an additional parameter b G {0, 1}. For matchings X(t, •), depending on b, 
this will either be equivalent to undoing the last step in the recursive definition of X (t, •) or 



yield (see Proposition |A.2) 



Definition A.l (Projection to t — 2 first vertices). Let t > 4 be an even integer. We define a 
function 

shrink t : (U 2 (U t ) U {0}) x {0, 1} -> U 2 (U t - 2 ) U {0} 

as follows. For b G {0, 1} we let shrinkt(0, 6) := 0. For M ^ 0, the definition is as follows; we 
let a(i) := ajiffi): 

1. {t - l,t - 2} e M: 

shrink (M, 1) := 

shrink t (M, 0) := M \ {{t - 1, t - 2}} 

2. {t - l,t - 3} e M: 

shrink t (M, 1) := (M \{{t-l,t- 3}, {t - 2, a(t - 2)}}) U {{t - 3, a(t - 2)}} 
shrink t (M,0) := 

3. {t - 2,t - 3} e M: 

shrink t (M, 6) := (M \ {{t - 1, a(t - 1)}, {t - 2,t - 3}}) U {{t - 3, a(t - 1)}} 

4. {t - 1, t - 2}, {t - 1, t - 3}, {t - 2, t - 3} g M: 

shrink t (Af, 1) := (M \ {{t - 1, a(t - 1)}, {t - 2, a(t - 2)}}) U {{a(t - 1), a(t - 2)}} 
shrink t (M, 0) := (M \ {{t - 1, a(t - 1)}, {t - 2, a(t - 2)}, {t - 3, a(* - 3)}}) U 
U {{t - 3, a(t - 2)}, {a(t - 1), a(t - 3)}} 

We omit the subscript t when it is clear from context. 

Proposition A. 2. The behavior of shrink^ for matchings X(t,ab) is a follows 



shrinkt(X(i, ab), c) 



X(t-2,a) ifb = c, 
else. 



Now we show how our shrink^ operation interacts with the fact whether or not two match- 
ings of Ut form a Hamiltonian cycle. Recall the n -operator: We have that the union of two 
perfect matchings M and M' of Ut is a Hamiltonian cycle if and only if M n M' = {Ut}, i.e., 
if the connectivity provided by the edges of the two matchings creates a single connected com- 
ponent that contains all vertices Ut- (We do not make further use of the meet operator but we 
find it convenient for specifying also the set of vertices on which we have a Hamiltonian cycle.) 

The main fact about this interaction is given by the following lemma. It shows that the fact 
whether perfect matchings M\ and M 2 form a Hamiltonian cycle depends directly on two pairs 
of projections of the two matchings. Note that the main statement of the lemma is given using 
addition modulo two. 
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Lemma A. 3. Let t > 4 be an even integer and let M±, M 2 £ Tl 2 (Ut). Let M? := shrink(Mj, j). 
Then 

[M x n M 2 = {U t }} = [Mi n Ml = {l/ t _ 2 }] + [M\ n M 2 ° = {tf t _ 2 }]. 

In other words, Mi U M 2 is a Hamiltonian cycle on Ut if and only if exactly one of M® U M 2 
and Ml U M$ is a Hamiltonian cycle on Ut- 2 . 

Proof. Unfortunately, there are quite a few cases that have to be considered. Each case is 
defined by the intersection of the set {{t — 1, t — 2}, {t — 1, t — 3}, {t — 2, t — 3}} with the two 
matchings; the cases are labeled by those of the edges that are contained in the two matchings. 
Due to symmetry via flipping the roles of Mi and M 2 a couple of cases are omitted. 

Let us clarify and recall a few things. The union of two perfect matchings on some vertex set 
is known to always give a disjoint union of cycles plus possibly a set of edges (the single edges 
are obtained when both matchings agree on the inclusion of some edges; we could regard them 
also as cycles of length two using a double edge). For the shr ink-operator it is only important 
which of the incident edges are in Mi and M 2 . For our subsequent arguments we also need to 
consider how the cycles in Mi and M 2 traverse the elements t — 1, — 2, and t — 3, and whether 
there are further connected components not containing any of these three elements. 

We let a := ckmj and (3 := o.m 2 be the corresponding functions for Mi and M 2 that for any 
element return the one that it is matched to. Accordingly, the cycles in Mi U M 2 will contain 
subpaths (a(t-l),t-l,(3(t-l)), (a(t-2), t-2, /3(t— 2)), and (a(*-3), i-3, /3(i-3)). Depending 
on the edges incident on t— 1, t— 2, and t — 3, we may get fewer paths: E.g., if {t — 1, t — 2} £ Mi 
then a(t - 1) = t-2 and a(t - 2) = t - 1; we get paths (/3(t - 2),i - 2,t - - 1) 

and (a(t-3),i-3,/3(t-3)). 

The shrink-operator will only change the edges with both endpoints in {t — l,t — 2,t — 
3, a(t — 1), a(t — 2), a(t — 3), /3(t — 1), /3(t — 2), f3(t — 3)}, i.e., all other edges are preserved in the 
projection/image. Thus, the way that the above paths are connected into cycles in Mi U M 2 is 
preserved in M® U M\ and M\ U M® ■ In particular, we may assume that Mi U M 2 contains no 
cycle disjoint from t — 1, t — 2, t — 3 since it would be preserved too; in that case neither of the 
three pairs can form a Hamiltonian cycle and there is nothing left to show. Otherwise, clearly, 
the way that a(t — 1), a(t — 2), ... , f3(t — 3) are connected affects directly whether Mi U M 2 is a 
Hamiltonian cycle. For the other two pairs of matchings created by shrink these connections 
are the same but the shrink-operator leads to different connections around t — 1, t — 2, and t — 3 
(in particular it removes t — 1 and t — 2), hence the total number of cycles may (and will) differ. 

1) {t — 1, t — 2} £ Mi n M2: Clearly Mi U M 2 is no Hamiltonian cycle as both matchings 
contain {£- l,i- 2}. We observe that Ml = shrink(Mi, 1) = and M\ = shrink(M 2 , 1) = 0. 
Hence neither M{* U M\ nor M\ U M® is a Hamiltonian cycle; this implies our claim. 

2) {i — 1, t — 3} £ Mi n -M2: Clearly Mi U M 2 is no Hamiltonian cycle as both matchings 
contain {t- l,t - 3}. We observe that Mf = shrink(Mi, 0) = and M$ = shrink(M 2 , 0) = 0. 
Hence neither M° U M\ nor U M® is a Hamiltonian cycle; this implies our claim. 

3) {t — 2, t — 3} £ .Mi n M2: Clearly Mi U M 2 is no Hamiltonian cycle as both matchings 
contain {t — 2, t = 3}. We will see that either both or none of M® U M 2 and M-J U M 2 form a 
Hamiltonian cycle. By Definition |A.l we have 



shrink t (M, b) := (M \ {{t - 1, a M (t - 1)}, {t - 2, t - 3}}) U {{t - 3, a M (t - 1)}}, 

when {i - 2,t - 3} G M. Thus Mj* = Mj and M\ = M 2 . It follows that M^ U M 2 X is a 
Hamiltonian cycle if and only if M\ U M 2 is a Hamiltonian cycle, as needed. 
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4) {t - l,t - 2} E Mi: We observe that Mi U M 2 contains subpaths (a(t- 3), i- 3, j3(t- 3)) 



and - l),t- l,t- 2,fi(t - 2)). By Definition |A.l| we get 
M? = M 1 \{{t-l,t-2}}, 

M\ = (M 2 \ {{t - 1, /3(t - 1)}, {t - 2, - 2)}}) U {{f3(t - 1), 0(t - 2)}}, 
Mi 1 = 0, 

Ml = (M 2 \ {{t - l,0(t - 1)}, {i - 2, 0(4 - 2)}, {t - 3, 0(t - 3)}}) U 
U {{t - 3, p(t - 2)}, - l),0(t - 3)}}. 

Hence, we know that M\ U M 2 cannot form a Hamiltonian cycle (and in the following, we 
will omit M| if its partner is the empty matching). Further, we get that M^ U M 2 contains 
subpaths (a(t — 3),i — 3,/3(i — 3)) and (J3(t — l),/3(i — 2)). Since these subpaths require the 
same configuration of the remaining vertices, as the subpaths of M\ U M 2 we get that M\ U M 2 
is a Hamiltonian cycle if and only if Mf U M 2 is a Hamiltonian cycle; this implies our claim. 

5) {i — 1, t — 3} G M%: We observe that M X UM 2 contains subpaths (a(t-2),t-2, f3(t-2)) 



and (/3(t — 1), £ — 1, t — 3,/3(t - 3)). By Definition |Al] we get 
M° = 0, 

Mf = (Mi \ {{* - 1, i - 3}, {t - 2, a{t - 2)}}) U {{t - 3, a(t - 2)}}, 
M 2 ° = (M 2 \ {{t - 1, P(t - 1)}, {t - 2, /3(t - 2)}, {t - 3, /3(t - 3)}}) U 
U {{t - 3, p(t - 2)}, {0(t - 1), /3(t - 3)}}. 

This time M{* U M 2 cannot form a Hamiltonian cycle, since M° = 0. For M-j 1 U M 2 we get 
subpaths (a(t — 2),t — 3,/3(t — 2)) and (f3(t — l),/3(t — 3)). Hence, the remaining edges com- 
plete Mi U M 2 to a Hamiltonian cycle if and only if the same is true for M^ U M$. 

6) {t - 2,t - 3} G Mi: We observe that M X UM 2 contains subpaths (a(t- 1), i- 1, /3(i- 1)) 



and (/?(£ — 2), £ — 2, £ — 3, /3(t — 3)). By Definition |A7L we get 



M° = (Mi \ {{t - 1, a(t - 1)}, {t-2,t- 3}}) U {{t - 3, a(t - 1)}}, 
M 2 = (M 2 \ {{t - 1, /3(t - 1)}, {t - 2, p{t - 2)}}) U {{/3(t - l),p(t - 2)}}, 
Ml = (Mi \ {{t - 1, a(i - 1)}, {t-2,t- 3}}) U {{t - 3, a(i - 1)}}, 
M 2 ° = (M 2 \ {{t - 1, /3(t - 1)}, {t - 2, /3(t - 2)}, {t - 3, /3(t - 3)}}) U 
U {{t - 3, /3(t - 2)}, {p{t - l)J(t - 3)}}. 

Now, for Ml U Ml we get subpaths (a(t -l),t- 3,/3(t - 3)) and (f3(t - l),(3(t - 2)), and 
for Ml U M 2 ° we get subpaths (a(t - 1), t - 3, /3(i - 2)) and (/3(t - l),j9(t - 3)). 

If the remaining edges contain any cycles, then none of Mi U M 2 , Mj 3 U M 2 , and M-j 1 U M 2 
form a Hamiltonian cycle (in the following we will tacitly ignore this case). Otherwise, all 
remaining edges form two vertex-disjoint paths that start and end in different vertices of {a(t — 
l),/3(t — l),f3(t — 2), f3(t — 3)}; for convenience, if, e.g., a(t — 1) = f3(t — 2) then we take this as 
a trivial path from a(t — 1) to j3{t — 2). There are three ways in which those paths can connect 
the four vertices: 

1. (a(t-l),...,/3(t-l)) and (/3(t - 2), . . . , (3{t- 3)): In this case Mi U M 2 is no Hamiltonian 
cycle, but both M® U M 2 and M\ U M 2 are Hamiltonian cycles. (This can be easily seen 
by combining the corresponding two subpaths with the two paths through the remaining 
edges: Either we get a single (Hamiltonian) cycle, or there are two cycles.) 
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2. (a(t -l),...,(3(t- 2)) and (/3(t - 1), . . . , /3(t - 3)): In this case M x U M 2 and M x ° U M\ 
are Hamiltonian cycles, and M\ U M 2 is no Hamiltonian cycle. 

3. (a(t -l),...,(3(t- 3)) and {j3(t - 1), . . . , fi(t - 2)): In this case M x U M 2 and U M 2 ° 
are Hamiltonian cycles, and M-f U M 2 is no Hamiltonian cycle. 

Thus, for all configurations of the remaining edges we get that M\ U M 2 forms a Hamiltonian 
cycle if and only if exactly one of U M 2 and M\ U M 2 is a Hamiltonian cycle. 

7) {t - 1, t - 2} £ Mi and {t - 1, t - 3} £ M 2 : We observe that Mi U M 2 contains a 



subpath (a(t -3),t-3,t- l,t- 2,(3(t - 2)). By Definition [ATI] we get 
M? = M 1 \{{t-l,t-2}}, 

Ml = (M 2 \{{ t -l,t- 3}, {t - 2, p(t - 2)}}) U {{t - 3, P(t - 2)}}, 
M\ = 0, 
M 2 ° = 0. 

Thus, M\ U M 2 is no Hamiltonian cycle. For M° U M 2 we get a single subpath (a(t — 3), t — 
3, — 2)). It follows easily that M\ U M 2 is a Hamiltonian cycle if and only if M® U M 2 is a 
Hamiltonian cycle. 

8) {t - 1, i - 2} G Mi and {£ — 2, t — 3} G M 2 : We observe that Mi U M 2 contains a 



subpath (a(t — 3), i — 3, t — 2, i — l,/3(t — 1)). By Definition A.l we get 



M 1 ° = M 1 \{{t-l,t-2}}, 

Ml = (M 2 \ {{t - 1, 0(t - 1)}, {t-2,t- 3}}) U {{t - 3, P(t - 1)}}, 
M] 1 = 0. 

Thus Ml U M 2 is no Hamiltonian cycle. For M^ U M 2 we get a single subpath (a(t — 3),i — 
3,/3(t - 1)). This implies that Mi U M 2 is a Hamiltonian cycle if and only if M-f U Mj is a 
Hamiltonian cycle. 

9) {t - 1, t - 3} £ Mi and {t — 2, t — 3} G M 2 : We observe that Mi U M 2 contains a 
subpath ( 



a/p/ia(t - 2),i - 2,t - 3,t - l,/3(t - 1)). By Definition [ATT] we get 
Mi = 0, 

Ml = (Mi \ {{* - 1, t - 3}, {t - 2, a(* - 2)}}) U {{t - 3, a(t - 2)}}, 
M 2 ° = (M 2 \ {{t - 1, 0(t - 1)}, {t - 2,t - 3}}) U {{t - 3, /3(i - 1)}}. 

Thus Mj 3 U M 2 is no Hamiltonian cycle. For Ml U M 2 we get a single subpath (a(t — 2),t — 
3, j3{t — 1)). This implies that Mi U M 2 is a Hamiltonian cycle if and only if Ml U M 2 is a 
Hamiltonian cycle. 

10) {{t - 1, t - 2}, {t - 1, t - 3}, {t - 2, t - 3}} n (Mi U M 2 ) = 0: In this final case 
we can only observe the trivial subpaths in Mi U M 2 , namely (a(t — 1), t — 1, /3(t — 1)), (a(t — 
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2),i-2,/3(i-2)), and (a(t - 3), t - 3, j3{t - 3)). By Definition [XT] we get 

M° = (Mi \ {{t - 1, a(t - l)},{t - 2, a(t - 2)},{t - 3, a(t - 3)}}) U 

U {{t - 3, a(t - 2)}, {a(t - 1), a(t - 3)}}, 
Ml = (M 2 \ {{t - 1, /3(t - 1)}, {i - 2, 0(i - 2)}}) U {{P(t - l),P(t - 2)}}, 
Ml = (Mi \ {{t - 1, a(t - 1)}, {t - 2, a(t - 2)}}) U {{a(t - 1), a(t - 2)}}, 
M 2 ° = (M 2 \ {{t - 1, 0(t - 1)}, {i - 2, /3(t - 2)}, {« - 3, /3(t - 3)}}) U 

U{{i-3,/3(i-2)},{/3(i-l),/3(i-3)}}. 

Thus, U M\ has subpaths 

(a(t - 1), a(t - 3)), (a(t - 2),t - 3, f3(t - 3)), (0(t - 1), f3(t - 2)), 
and Mi 1 U M^ has subpaths 

(a(t - 1), a(t - 2)), {a{t - 3), t - 3, f3(t - 2)), (/3(t - 1), f3(t - 3)). 

Thus, whether any of the three unions of perfect matchings is a Hamiltonian cycle depends 
again on the paths given by all remaining edges. As before, if those edges give rise to cycles 
disjoint from t — 1, t — 2, t — 3 then none of the three unions is a Hamiltonian cycle and we are 
done. Else, the remaining edges serve exactly to create three vertex disjoint paths between the 
vertices a(t — 1), a(t — 2), a(t — 3), f3(t — 1), j3{t — 2), and /3(t — 3); there are 15 configurations of 
such paths (i.e., partitions of this set of endpoints into three pairs of connected vertices). Note 
that again we treat cases like a(t — 1) = /3(t — 2) has a path connecting a(t — 1) and /3(t — 2). 

It can be easily verified that for each of the 15 configurations we get the desired property 
that Mi U M 2 is a Hamiltonian cycle if and only if exactly one of M® U M\ and M\ U M® is a 
Hamiltonian cycle. For completeness, we provide a table containing all 15 cases; see Table [T] 
This completes the proof. □ 



As a special case of Lemma A. 3 when the second matching is in our intended basis, we get 
the following lemma. Both lemmas are used for proving the main theorem of this section. 

Lemma A. 4. Let t be an even integer t > 6 and let M £ n 2 (Ut). Furthermore let X = X(t, ab) 
where b £ {0, 1} and a denotes a 0, 1-string of length | — 2. Then MUX is a Hamiltonian cycle 
if and only if shrink(M, b) U X(t — 2, a) is a Hamiltonian cycle, i.e., 

[M n X(t, ab) = {U t }} = [shrink(M, b) n X(t - 2, a) = {U t - 2 }}. 



We are now set up to prove Theorem 3.4 let us recall the theorem statement first. 



Theorem A. 5. 3.4 Let t > 2 be an even integer and let Ml,M 2 £ n 2 (C/t). It holds that 
[Mi n M 2 = {u t }] = t Ml n x ^ a ) = Ml • t M2 n x ^ s ) = t 17 *}]' 

aejO,!}'/ 2 - 1 



where X(t,a),X(t,a) £ X t according to Definition 3.1. (Each matching in X f occurs exactly 
twice, once as X(t,a) and once as X(t,a).) 

We establish a factorization of the Matching Connectivity matrix as the product of two 
smaller rectangular matrices: The first matrix has rows labeled by (all) perfect matchings and 
columns labeled by basis matchings X(t, a). The second matrix has rows labeled basis match- 
ings X(t,a) and columns labeled by (all) perfect matchings. The entries of both rectangular 
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config 


juration of 


remaining edges 


as paths 


between 




Mi U M 2 


AfJ U 


Ml U M 2 ° 


a(t — 


1), a{t - 2 


), a(i - 


-3),0(t- 


-i),p(t- 


-2), p{t 


-3) 


is H-cycle 


is H-cycle 


is H-cycle 


a(t - 
a(t - 
a(t- 


l)-a(t-2) 
l)-a(t-2) 
l)-a(t - 2) 


, a(t — 
, a(t — 
, a(t — 


3)-/3(t - 
3)-/9(t - 
3)-/9(* - 


!),£(*- 

2) ,/3(t- 

3) , /3(t - 


2)-/3(i- 


3) 
3) 
2) 


1 
1 




1 
1 








a[t — 


l)-a(t - 3) 


, a(t — 


2)-/3(* - 


l),/9(t- 


2)-P(t - 


3) 


1 





1 


a(t — 


l)-a(t-3) 


, a(t — 


2)-/3(i - 


2),/3(t- 


1)-P(t - 


3) 











a(t — 


l)-a(t - 3) 


, a(t — 


2)-/3(i - 


3), - 


l)-/3(t - 


2) 


1 





1 


a(t - 
a(t — 
a(t — 


l)-(3(t - 1) 
l)-/3(t-l) 
l)-/3(i-l) 


, a(t — 
, a(t — 
, a[t — 


2)-a(t - 
2)-p(t - 
2)-0(t - 


3), (3(t- 

2) , a(t - 

3) , a(t - 


2) -(3(t - 

3) -P(t - 
3)-(3(t - 


3) 
3) 
2) 







1 
1 



1 
1 




a(t — 
a(t — 
a(t - 


!)-/?(* -2) 
- 2) 
l)-0(t - 2) 


, a[t — 
, a(t — 
, a(t — 


2)-a(t - 
2)-/3(t - 
2)-/3(t - 


3),(3(t- 
1), a(t - 
3), a(t - 


l)-/3(t - 
3)-/3{t - 
3)-(3(t - 


3) 
3) 
1) 


1 



1 


1 
1 





1 
1 


a(t - 
a(t - 
a(t — 


- 3) 
l)-/3(i - 3) 
l)-/5(t - 3) 


, a(t — 
, a(t — 
, a[t — 


2)-a(t - 
2)-p(t - 
2)-P{t - 


3),(3(t- 

1) , a(t - 

2) , a(t - 


l)-(3(t - 
3)-/3(t - 
3)-(3(t - 


2) 
2) 
1) 


1 
1 





1 
1 


1 



1 



Table 1: All 15 configurations for paths connecting the vertices a(t — 1), a(t — 2), a(t — 3), (3(t — 
1), f3(t - 2), and (5{t - 3), and their effect on whether Mi U M 2 , Mf U M 2 \ and M x x U M 2 ° form 
a Hamiltonian cycle (denoted by 1 and otherwise. 



matrices are as in the large one: we have 1 if the two matchings form a Hamiltonian cycle and 
else (so in fact they are submatrices of the large matrix). Note that the two perfect matchings 
from Xt in each summand are exactly the unique pairs that give Hamiltonian cycles (among 
basis matchings). Also note that the main statement of the theorem is given in arithmetic 
modulo two. 

Proof. Let us quickly verify the theorem for the base case t = 2. There is a unique perfect 
matching {{0, 1}} for Ui = {0, 1}. Thus, we can only pick M\ = M 2 = {{0, 1}}. Hence, for the 
left hand side of the theorem statement we get 

LHS = [Mi n M 2 = {U 2 }\ = 1. 

Regarding the right hand side the only choice for a G {0, 1}*/ 2_1 in the sum is a = e (of 
course e = e). We recall that X(2, e) = {{0, 1}}, thus we get 

RHS = [M x n X(2,e) = {U 2 }] ■ [M 2 n X(2,e) = {U 2 }\ = 1. 

This completes the base case t = 2. 

Now, we will perform an inductive argument to show that the theorem holds also for all 



even t > 4. We start from the right hand side of the theorem statement and apply Lemmas A. 3 



and A. 4 to transform it into the left hand side. As a first step, we split the right hand side 



sum into two parts, depending on the last position of a (to do so formally we shorten a by one 
position and explicitely specify or 1 to be appended). 

RHS = Yl [ Ml n a0 ) = " t M 2 n a0) = i U t}] 

ae{0,l} t/2 ~ 2 

+ Yl [M 1 nx{t,ai) = {u t }]-[M 2 nx{tM) = {u t }} 

aSfO,!}'/ 2 " 2 
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Note that aO = al and al = aO. Let us recall the statement of Lemma A. 4 



[M n X(t,ab) = {U t }] = [shrink t (M,£>) n X(t -2,a) = {C/ t _ 2 }]. 
Using the shorthand M- := shrinkj(Mj, j) we can thus rewrite the RHS as follows. 

RHS = i M i n X ( l ~ 2 ' a ) = {^-2}] • i M 2 n X(t - 2, a) = {U t - 2 }] 

ae{0,l}* /2 " 2 

+ Yl [Mi 1 ni(f-2,o) = {M'Wn^-2,«) = H 

a6{0,l} i/2 - 2 

We note that both sums go over all bitstrings a of length ^ — 2 = — 1. Thus, after brief 
inspection of the summands, we see that we can apply the inductive assumption to both sums 
since m{ £ n 2 (£/t_ 2 ). We obtain 

RHS = [M° n Ml = {U t - 2 }] + [Ml n M 2 ° = {U t - 2 }}. 



Finally, by a direct application of Lemma A. 3 we obtain 

RHS = [Mi n M 2 = {U t }} = LHS, 

which is what we intended to prove. Thus, the theorem holds for all even integer t > 2, as 
claimed. □ 
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